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Abstract 

We formulate a complete theory of Edge Radiation based on a novel method 

relying on Fourier Optics techniques. Similar types of radiation like Transition 
Undulator Radiation are addressed in the framework of the same formalism. Special 
attention is payed in discussing the validity of approximations upon which the 
theory is built. Our study makes consistent use of both similarity techniques and 
comparisons with numerical results from simulation. We discuss both near and 
far zone. Physical understanding of many asymptotes is discussed. Based on the 
solution of the field equation with a tensor Green's function technique, we also 
discuss an analytical model to describe the presence of a vacuum chamber. In 
particular, explicit calculations for a circular vacuum chamber are reported. Finally, 
we consider the use of Edge Radiation as a tool for electron beam diagnostics. 
We discuss Coherent Edge Radiation, Extraction of Edge Radiation by a mirror, 
and other issues becoming important at high electron energy and long radiation 
wavelength. Based on this work we also study the impact of Edge Radiation on 
XFEL setups and we discuss recent results. 
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1 Introduction 



Synchrotron Radiation (SR) sources from bending magnets are brilliant, and 
cover the continuous spectral range from microwaves to X-rays. However, 
in order to optimally meet the needs of basic research with SR, it is desir- 
able to provide specific radiation characteristics, which cannot be obtained 
from bending magnets, but require special magnetic setups, called inser- 
tion devices. These are installed along the particle beam path between two 
bending magnets, and introduce no net beam deflection. Therefore, they can 
be incorporated in a given beamline without changing its geometry. Undu- 
lators are a typical example of such devices, generating specific radiation 
characteristics in the short wavelength range. 

The history of SR utilization in the long wavelength region (from microm- 
eter to millimeter) is more recent than that in the short wavelength range. 
Long wavelength SR sources may have a strong potential for infrared spec- 
troscopy or imaging techniques. In fact, they are some order of magnitude 
brighter than standard thermal sources in the same spectral range. 

Large angles are required to extract long wavelength SR from bending 
magnets, because the "natural" opening angle in this case increases up to 
several tens milliradians in the far-infrared range. However, the situation 
changes dramatically if a straight section is introduced between two bends, 
like in Fig. 1(a). Long- wavelength radiation emitted by relativistic electrons 
in this setup is called Edge Radiation (ER), and presents a significantly 
smaller opening angle than standard SR from bends (see, among others, [1]- 
ifl^ ). In other words, in the long wavelength region (compared to the critical 
bending-magnet radiation wavelength) a simple straight section between 
bends can play the role of a kind of insertion device. 

ER and bending magnet radiation have equivalent brightness. In fact, the 
physical process of ER emission is not different from that of radiation emis- 
sion from a single bend. However, radiation from the setup in Fig. 1(a) 
exhibits special features. Due to a narrower opening angle of ER over SR 
from bends, as the wavelength gets longer ER yields significant advantages 
in terms of simplicity of the photon beamline. 

ER theory is a part of the more general SR theory, very much like Undulator 
Radiation (UR) theory is a part of SR theory. Similarly to the UR case, also for 
ER the knowledge of the applicability region of the far-field formulas and 
corrections for near-field effects are of practical importance. In most practical 
cases, the distance between ER source and observer (i.e. the first optical 
element of the photon beamline) are comparable or even much smaller than 
the length of the straight section, which plays the role of the length of the 
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Fig. 1. Four main types of edge radiation setups: (a) Far-infrared beamline for 
synchrotron radiation source using edge radiation, (b) Arrival-time monitor for 
XFEL source using optical coherent edge radiation, (c) Electron bunch length mon- 
itor for XFEL using far-infrared coherent edge radiation, (d) Ultra-short electron 
bunch diagnostic for laser-plasma accelerator facility using optical coherent edge 
radiation. 
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insertion device for ER. 



In this paper we developed a theory of near-field ER based on Fourier Optics 
(FO) techniques. These techniques can be exploited without limitations for 
ER setups, because the paraxial approximation can always be applied in 
the case of electrons in ultra-relativistic motion ||l5(l . The use of the paraxial 
approximation allows reconstruction of the field in the near-zone from the 
knowledge of the far-field data. The solvability of the inverse problem for the 
field allows characterization of any ER setup, starting from the far-zone field, 
in terms of virtual sources. These sources exhibit a plane wavefront, and can 
be pictured as waists of laser-like beams. Using this kind of description we 
develop our theory in close relation with laser-beam optics. In particular, 
usual FO can be exploited to characterize the field at any distance, providing 
a tool for designing and analyzing ER setups. 

It is the purpose of this article to discuss the principles of production and 
properties for all applications of ER. First, we treat the relatively simple case 
of ER from a setup composed by straight section and two bending magnets 
at its ends (see Fig. 1(a)). We begin calculating an analytical expression for ER 
from a single electron in the far-zone. Then, we characterize the near-zone 
with the help of the virtual-source technique. Two alternative techniques 
for the field propagation are given, based on a single virtual source located 
in the middle of the ER setup, and based on two virtual sources located at 
its edges. 



Based on this study-case we turn to analyze a more complicated setup, 
consisting of an undulator preceded and followed by two straight sections 
and two bends (see Fig. 1(b)). ER from this kind of setup is commonly known 
as Transition Undulator Radiation (TUR). The first study on TUR appeared 
more than a decade ago in |fl6t]. In that work it was pointed out for the first 
time that, since an electron entering or leaving an undulator experiences 
a sudden change in longitudinal velocity, highly coUimated radiation with 
broadband spectrum, similar to transition radiation, had to be expected 
in the low-frequency region in addition to the usual UR. Reference | |l6|l 
constitutes a theoretical basis for many other studies. Here we remind only 
a few ^ [l^ 18, 19, |2^, dealing both with theoretical and experimental 
issues. More recently, TUR has been given consideration in the framework 
of X-ray Free-Electron Laser (XFEL) projects like 121, 22, 23l]. For example, 
an arrival-time monitor for XFELs using infrared coherent ER from a setup 
similar to that in Fig. 1(b) has been proposed in [24], which should be used 
for pump-probe experiments with femtosecond-scale resolution. In view of 
these applications, there is a need to extend the characterization of TUR to 
the near-zone, and to the coherent case. From this viewpoint, specification 
of what precedes and follows the undulator is of fundamental importance. 
As has been recognized for TUR many years ago |i90, if this information 
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is not known, any discussion about the intensity distribution of TUR is 
meaningless. According to our approach, the two straight sections and the 
undulator in the setup in Fig. 1(b) will be associated to virtual sources with 
plane wavefronts. The field from the setup can then be described, in the near 
as well as in the far-zone, as a superposition of laser-like beams, radiating 
at the same wavelength and separated by different phase shifts. 

Our study makes consistent use of both dimensional analysis and com- 
parisons with outcomes from numerical simulation. All simulations in this 
paper are performed with the help of the computer code SRW [25] . 

There are, however, situations when existing computer codes cannot predict 
the radiation characteristics. One of these is the case when perturbations of 
the long-wavelength radiation by vacuum chambers are present, which may 
potentially affect the performance of ER setups. Since the diffraction size of 
the THz radiation exceeds the vacuum chamber dimensions, characteriza- 
tion of far-infrared ER must be performed accounting for the presence of a 
waveguide. In order to deal with this situation we developed a theory of ER 
in a waveguide. The task to be solved differs from the unbounded-space case 
only in the formulation of boundary conditions. The paraxial approxima- 
tion applies as in the unbounded-space case. Only, on perfectly conducting 
walls the electric field must be orthogonal to the vacuum chamber surface. 
As in the unbounded-space case, one can use the Green's function approach 
to solve the field equations. The presence of different boundary conditions 
complicates the solution of the paraxial equation for the field, which can 
nevertheless be explicitly found by accounting for the tensorial nature of 
the Green's function. Here we take advantage of a mode expansion approach 
to calculate ER emission in the metallic waveguide structure. We solve the 
field equations with a tensor Green's function technique, and we extract 
figure of merits describing in a simple way the influence of the vacuum 
chamber on the radiation pulse as a function of the problem parameters. We 
put particular emphasis on a vacuum chamber with circular cross-section, 
which is natural for future linac-based sources (XFELs and Energy Recovery 
Linacs (ERLs) 126J). 

Finally, we address the long-standing interest of electron beam characteriza- 
tion for linac-based sources and laser-plasma accelerators. One possibility 
to perform electron-beam diagnostics in these kind of facilities is to use 
coherent ER. This is an attractive tool, because it can provide valuable and 
detailed information on the electron beam. By detecting coherent ER, 3D dis- 
tributions, divergence, micro-bunching may be measured in principle (See 
Figs. 1(c) and (d)). Usually, electrons in accelerators are highly collimated 
and monochromatic. In this case, coherent ER can be used for longitudinal 
and transverse beam-size monitoring. In contrast to this, electrons gener- 
ated in laser-plasma interactions have different properties compared with 
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those in conventional accelerators. Namely, in this case, electrons have both 
divergence angle and energy distribution. We address applications of coher- 
ent ER by studying, first, the relatively simple case when only the influence 
of longitudinal and transverse structure factor of the electron bunch is ac- 
counted for. In particular, we give particular attention to the case when 
microbunching at optical wavelengths is imprinted onto an electron bunch, 
and we analyze the more complicated case of a bunch produced in a laser- 
plasma accelerator, i.e. accounting for divergence angle distribution of the 
beam. Energy distribution can be easily accounted for, based on this analy- 
sis. The problem of extraction of ER by a mirror, strictly related to diagnostics 
applications, is included too. 



2 General relations for edge radiation phenomena 

2.1 Physical discussion of some numerical experiment 

This Section constitutes an attempt to introduce ER theory to readers in 
as intuitive and simple a fashion as possible by simulating the angular 
spectral flux as a function of observation angles for the geometry in 1(a). 
For this purpose we take advantage of the code SRW | | 25|] , which provides 
a numerical solution of Maxwell's equations. 

The origin of a Cartesian coordinate system is placed at the center of straight 
section. The z-axis is in the direction of straight section and electron motion 
is in xz plane. Parameters of the problem are the radiation wavelength A, the 
radius of the bend K, the relativistic Lorentz factor y, the length of the straight 
section L and, additionally, the position of the observation plane down the 
beamltne, z. We work in the far zone. In this Section it is operatively defined 
as a region where z is large enough, so that the simulated angular spectral 
flux does not show dependence on z anymore. 

ER carries advantages over bending magnet radiation in the limit for A/ » 
1, where Ac ~ -R/y'^ (here A = A/(27i) is the reduced wavelength) is the 
critical wavelength of bending magnet radiation. We will work, therefore, 
in this limit. We set y = 3.42 • 10"^ (17.5 GeV), R = 400 m, which are typical 
values for XFELs. Note that in this case Ac - O.IA. Here we take A = 400 
nm. We begin with the case L = (bending magnet), and we increase the 
straight section length (see Fig. 2). As one can see from the figure, radiation 
becomes more collimated, up to about L - y^A ^ 100 m (case (d)), where 
the collimation angle reaches 1/y ~ 30/jrad. Further increase of L only leads 
to the appearance of finer structures in the radiation profile. It is important 
to remark that the total number of photons in the ±1 mrad window shown 
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in Fig. 2(a) is roughly the same in the ±100jUrad window in Fig. 2(d). It is 
clear that the length of the straight section L is strongly related with the 
collimation of the radiation. 



2.2 Similarity techniques 



To study ER further we apply similarity techniques. Similarity is a special 
symmetry where a change in scale of independent variables can be com- 
pensated by a similarity transformation of other variables. This is a familiar 
concept in hydrodynamics, where the cardinal example is given by the 
Reynolds number. Similarity allows one to reduce the number of parame- 
ters to a few dimensionless ones that are directly linked to the physics of the 
process, and that control it in full. Such parameters are found by analysis 
of the underlying equations characterizing the system under study. In this 
Section we limit ourselves to list them, to show their correctness with the 
help of the code SRW, and to describe their physical meaning. This allows 
one to obtain general properties of the ER process. A comprehensive theory 
of ER will be presented in the following Sections. 

For the setup in Fig. 1(a), two dimensionless parameters controlling the 
radiation characteristics can be extracted from Maxwell's equations. In the 
next Section we will show how these parameters can be derived. Here we 
limit ourselves to write them: 



^ = — ' "^ = 7^- 

The detector is supposed to be far away from the source so that the above- 
given definition of far-zone holds. 

The most important general statement concerning ER is that all possible sit- 
uations correspond to different values of the two dimensionless parameters 
6 and (p. 

Note that the working limit A /Ac » 1 means (/) • 6 «; 1 in terms of dimen- 
sionless parameters. For any two cases characterized by the same values of 
6 and (p the angular spectral flux from set up in Fig. 1(a) will 'Took" the same 
in terms of angles scaled to VA/L, i.e. 9 = 9/ s/A/L. In other words, data 
for different sets of problem parameters corresponding to the same values 
of 6 and (p reduce to a single curve when properly normalized. We tested 
the scaling properties of ER by running numerical simulations with the first 
principle computer code SRW. We used two different sets of dimensional 
parameters corresponding to the same case in terms of parameters 6 and cp, 
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Fig. 2. Illustrative calculations of the effect of bending magnet separation on the 
directivity diagram of the radiation. The bending magnet radius R = 400 m, the 
relativistic factor y - 3.42 • 10^, and the wavelength of interested A = 400 nm are 
fixed, while the straight section length varies from L = up to L » y^A ^ 100 
m. In this setup (as well as in all others in this paper) A » Ac - 0.1 A. Case (a) is 
a bending magnet setup. Case (b) is a complex setup, where the radiation beam 
divergence is practically the same as in (a). Case (c) illustrates an ER setup. Bending 
magnet separation dramatically lowers the radiation beam divergence, (d) Optimal 
bending magnet separation. The straight section length L ^ y^A corresponds to a 
radiation beam divergence G ^ 1/y. (e) Further increase of L only leads to the 
appearance of finer structures in the radiation profile. 2D plots on the left show 
the angular spectral flux as a function of the horizontal and vertical angles Gx and 
Gy for various lengths of the straight section. Middle plots are obtained cutting 
the 2D angular spectral flux profile at x = 0. Right plots show a schematic of the 
considered layout. 
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Fig. 3. Verification of similarity techniques. Left and right plots show the normalized 
angular spectral flux as a function of the horizontal and vertical angles Ox and 6y 
respectively (at = and 0^ = respectively), (a) Case 6 ^ 0.43 and cp ^ 6.7 ■ 
Solid curve is the result of SRW calculations with L - Q.5 m, R - 400 m, A = 400 
nm at 17.5 GeV. Dotted curve is the result for L = 1 m, R = 800 m, A = 800 nm at 
17.5 GeV. (b) Case 6 «; 1 and c/) ^ 4. Solid curve is the result of SRW calculations 
with L = 300 m, R = 400 m, A = 400 nm, at 17.5 GeV (corresponding to 6 ^ 7 • lO"'^). 
Dotted curve is the result for L = 150 m, R = 400 m, A = 800 nm at 8.5 GeV 
(corresponding to 6 ^ 2 • 10"^). 




Fig. 4. Illustration of self-similarity techniques. Left and right plots show the nor- 
malized angular spectral flux as a function of the horizontal and vertical angles Ox 
and 0y respectively (at Oy - ^ and 0x = respectively). The angular spectral flux 
profile asymptotically approaches the self-similar form I/Imax = F{Ox, Oy) for 6 «c 1 
and (p <^\. Solid curve is the result of SRW calculations with 6 ^ 0.02 and (p ^ 0.13. 
Dotted curve is refers to the case 6 = 0.01 and (p - 0.27 instead. 

and we checked that the angular spectral fluxes normalized to their max- 
imal values are identical. Results are presented in Fig. 3 and Fig. 4, where 
the normalized angular spectral flux is indicated with I/Imax- 

When 6 ~ 1, the presence of the bending magnet radiation strongly influ- 
ences the radiation profile (see Fig. 3(a)). When 6 decreases up to 6 <c 1, one 
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can neglect bending magnet contributions (see Fig. 3(b)): what is left in this 
case is ER. These situations are realized, for example, if one works at fixed A, 
y and R while increasing the length L as in the case of Fig. 2. It follows that 
6 is responsible for the relative weight of ER and bending magnet radiation 
contributions in the radiation profile. Since we are interested in ER emis- 
sion, it is natural to consider more in detail the limit for 6 «; 1. In this case, 
results are independent on the actual value of 6, and the only parameter left 
is (p. This fact can be seen from Fig. 3(b), where the two sets of dimensional 
parameters refer to two different value of 5 <c 1. We will name this situation 
the sharp-edge asymptote. 

In the limit for « 1, the opening angle of the radiation is independent 
of the actual value of (p too. In this case we obtain the universal plot 
shown in Fig. 4, and one talks about a self-similar behavior of the angu- 
lar spectral flux profile, which asymptotically approaches the self-similar 
form ///max = F{dx, dy)- Notc that the separation distance L between the 
bends dramatically lowers the radiation beam divergence, but the charac- 
teristic angle of emission is still larger than 1 /y. In fact, radiation peaks at 
6 - 2.2 VA/L. When (p increases, radiation becomes better and better coUi- 
mated, up to angles ~ l/y. This happens for values (p - 1. Radiation has 
reached the best possible collimation angle and further increase of (p (see 
Fig. 3(b)) only modifies fine structures in the radiation profile. 



2.3 Qualitative description 



It is possible to present intuitive arguments to explain why all problem 
parameters (J^, y, L and A) are effectively grouped in 6 and (p. 

To this purpose let us consider first the parameter 6. By definition, 1/6 is a 
measure of the straight section length L in units of a characteristic length 



To explain the meaning of the quantity '\IR^A, following \Zn\ we consider 
Fig. 5(a), and we focus on the region of parameters A <^ R and » 1. 
A posteriori, this region of parameters will turn out to correspond to an 
angular dimension along the trajectory 20 <?c 1 within the bending magnet. 
Radiation from an electron passing through the setup is observed through 
a spectral filter by a fixed observer positioned on the tangent to the bend 
at point P. Electromagnetic sources propagate through the system, as a 
function of time, as shown in Fig. 5(b). However, electromagnetic signals 
emitted at time t' at a given position x{t') arrives at the observer position 
at a different time t, due to finite speed of light. As a result, the observer in 
Fig. 5(a) sees the electromagnetic source motion as a function of t. What one 
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Fig. 5. Geometry for SR from a bending magnet. 

needs to know, in order to calculate the electric field, is the apparent motion 
x{t) shown in Fig. 5(c), which is a hypocycloid, and not the real motion 
x{t'). In fact, the electric field at the observation point is proportional to the 
second derivative of the x-coordinate with respect to the retarded time f, 
because the observer sees everything as delayed. We discuss the case when 
the source is heading towards the observer. Using the fact that <^c 1, one 
obtains the well-known relation dt/dt' = 1/2 • (1/y^ + 6^). The observer sees 
a time-compressed motion of the sources, which go from point A to point 
B in an apparent time corresponding to an apparent distance lRddt/{dt'). 
Let us assume (this assumption will be justified in a moment) 6^ > In 
this case one has 2R6dt/{dt') - Rd^. Obviously one can distinguish between 
radiation emitted at point A and radiation emitted at point B only when 
R6^ » A, i.e. for 6 » (A/Ry^^. This means that, as concerns the radiative 
process, we cannot distinguish between point A and B on the bend such 
that Rd < {R'^Ay^^. It does not make sense at all to talk about the position 
where electromagnetic signals are emitted within Lfi, = (R^A)^^^ (here we 
are assuming that the bend is longer than Lfi,). This characteristic length 
is called the formation length for the bend. The formation length can also 
be considered as a longitudinal size of a single-electron source. Note that 
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a single electron always produces diffraction-limited radiation • A0 ~ A, 
d being the transverse size and A6 the divergence of the source. Since d ~ 
Lfi,Ad, it follows that the divergence angle A6 is strictly related to Lfh and A: 
e ~ 7X717. One may check that, using Lfi, ~ yfWX^, one obtains 6 ~ ^XjR; 
in particular, at A ~ Ac ~ R/y^ one obtains 6 ~ I/7, as is well-known for 
bending magnet radiation. 

Let us now consider the case of a straight section of length L inserted between 
the two halves of a bend. Since we cannot distinguish between points within 
Lf},, the case L = is obviously indistinguishable from the case L <^ Lf. 
Significant deviations from the bending magnet case are to be expected 
when L > Lfj,, i.e. when 6 < 1. This hints to the fact that 6 is responsible 
for the relative weight of ER and bending magnet radiation contributions 
in the radiation profile. 

Let us now discuss the parameter (p. By definition, (p is a measure of the 
straight section length L in units of a characteristic length y^A. One can still 
use the same reasoning considered for the bend to define a region of the 
trajectory where it does not make sense to distinguish between different 
points. In the case of a straight section of length L connecting A and B, 
dt/dt' = l/{2y^). It follows that the apparent distance AB is equal to L/(2y^). 
Since it does not make sense to distinguish between points within the ap- 
parent electron trajectory such that L/{2y^) < A, one obtains a critical length 
of interest ~ y^A. This hints to the fact that for values (p ^ 1 radiation has 
reached the best collimation angle. 

Note that for ultrarelativistic systems in general, the formation length is al- 
ways much longer than the radiation wavelength. This is related with a large 
compression factor dt/dt'. For comparison, in the case of non-relativistic mo- 
tion the compression factor dt/dt' - 1, and the formation length is simply of 
order of the radiation wavelength. The counterintuitive result follows, that 
for ultrarelativistic systems one cannot localize sources of radiation within 
a macroscopic part of the trajectory. 



3 Paraxial approximation 

In the next two Sections we present a complete theory of ER. All electrody- 
namical theories are based on the presence of small or large parameters. 

In general, the theory of Synchrotron Radiation (SR) is based on the ex- 
ploitation, for ultra-relativistic particles, of the small parameter y~^. By this. 
Maxwell's equations are reduced to much simpler equations with the help 
of paraxial approximation. ER theory constitutes a particular case of SR 
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theory, based on the extra small-parameter 6. 

Here and everywhere else in this paper we will make consistent use of 
Gaussian units. 

In this Section we deal with the paraxial approximation of Maxwell's equa- 
tions. We will treat both near and far zone cases, with special attention to 
the applicability region of equations describing ER in different regions of 
the parameter space. 

Whatever the method used to present results, one needs to solve Maxwell's 
equations in unbounded space. We introduce a cartesian coordinate sys- 
tem, where a point in space is identified by a longitudinal coordinate z and 
transverse position r. Accounting for electromagnetic sources, i.e. in a re- 
gion of space where current and charge densities are present, the following 
equation for the field in the space-frequency domain holds in all generality: 

c^V^I + (o^E = 4:nc^Vp - 4:nia)J , (2) 



where p(z, r, co) and j{z, r, oS) are the Fourier transform^ of the charge den- 

sity p(z, r, t) and of the current density /(z, r, t). Eq. (2) is the well-known 

Helmholtz equation. Here E indicates the Fourier transform of the electric 
field in the space-time domain. 

A system of electromagnetic sources in the space-time can be conveniently 
described by p{z,r,t) and j{z,r,t). Considering a single electron and using 
the Dirac delta distribution, we can write 

p (z, r*, t) = -e6 (r*- ^oit)) 6{z - Zo(0) = (r - ro(z)) dl'-^-t] (3) 

f{z,r,t) = x^it)p{z,r,t) , (4) 



where (zo(0/j^o(0) ^rid v{t) are, respectively, position and velocity of the 
particle at a given time Hn a fixed reference frame, Vz is the longitudinal 
velocity of the electron , and {-e) is the electron charge. Additionally, we 

^ We explicitly write the definitions of the Fourier transform and inverse transform 
of a function f{t) in agreement with the notations used in this paper. The Fourier 
transform and inverse transform pair reads: 

f{co) = I dt fit) exp [icot] ; f{t) = I dco /(co) exp [-icot] . 
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defined the curvilinear abscissa s(z) - vt{z), where v = \v{t{z))\ is a constant. 
In the space-frequency domain the electromagnetic sources transform to: 

p{r, z, co) = —6 (r- ro(z)) exp 

"VzyZ) 



icos{z) 



V 



(5) 



and 



j{r,z,a)) = v{z)p{f,z,co) 



(6) 



Since we will only be interested in the transverse components of the field, 

-» 

from now on we will consider the transverse field envelope E, a 2D vector 

defined in the space-frequency domain as £ = E_iexp[-icoz/c], the sym- 
bol "_L" indicating projection on the transverse plane. By substitution in 
Heknholtz equation we obtain 



+ 



2ia) d 



E = —— exp 



ICO 



'siz) 



V 



C \ f7 



6(r'-Uz))i7) 



where, according to our notation, Vx indicates a gradient with respect to 
transverse coordinates only, and v± is the transverse velocity of the electron. 
Eq. (7) is still fully general and may be solved in any fixed reference system 
{x, y, z) of choice with the help of an appropriate Green's function. 

When the longitudinal velocity of the electron, Vz, is close to the speed of 
light c, one has yl :^ 1, where = 1 - vl/c^. The Fourier components 
of the source are then almost synchronized with the electromagnetic wave 
travelling at the speed of light. Note that this synchronization is the reason 
for the time compression factor described in Section 2.3. In this case, the 
phase co{s{z)/v - z/c) is a slow function of z compared to the wavelength. 
For example, in the particular case of motion on a straight section, one has 
s(z) = z/Vz, so that co{s{z)/v - z/c) = coz/{2ylc), and if » 1 such phase 
grows slowly in z with respect to the wavelength. For a more generic motion, 
one similarly obtains: 



CO 



's{z2) - s(zi) 




2yl{z)c ' 



(8) 



Mathematically, the phase in Eq. (8) enters in the Green's function solution 
of Eq. (7) as a factor in the integrand. As we integrate along z, the factor 
co{s{z)lv-zlc) leads to an oscillatory behavior of the integrand over a certain 
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integration range in z. Such range can be identified with the value of Z2 - Z\ 
for which the right hand side of Eq. (8) is of order unity, and it is naturally 
defined as the radiation formation length of the system at frequency 
(jj. Of course there exist some freedom in the choice of such definition: 
"order of unity" is not a precise number, and reflects the fact that there is 
no abrupt threshold between "oscillatory" and "non-oscillatory" behavior 
of the integrand in the solution of Eq. (7). It is easy to see by inspection of 
Eq. (8) that if is sensibly smaller than c, but still of order c, i.e. ~ c 
but l/y^ ~ 1/ then Lf ~ K. On the contrary, when t^z is very close to c, i.e. 
\ly\<^\, the right hand side of Eq. (8) is of order unity for Ly^ = Z2 - Zi » A. 

When the radiation formation length is much longer than A, £ does not vary 
much along z on the scale of A, that is | <9z£x,y cf/c | Ex,y |. Therefore, the 
second order derivative with respect to z in the operator on the left hand 
side of Eq. (7) is negligible with respect to the first order derivative, and Eq. 
(7) can be simplified to 



+ 



lioj d \ ^ 4ne 



c dzi 



exp 



10) \ 



's(z) 



-7U±(z) - Vj 



5 (r- fo{z)) , 



(9) 



where, as said before, we consider transverse components of E, and we 
substituted Vz{z) with c, based on the fact that 1/yl <^ 1. Eq. (9) is Maxwell's 
equation in paraxial approximation. Eq. (7), which is an elliptic partial dif- 
ferential equation, has thus been transformed into Eq. (9), that is of parabolic 
type. Note that the applicability of the paraxial approximation depends on 
the ultra-relativistic assumption » 1 but not on the choice of the z axis. 
If, for a certain choice of the longitudinal z direction, part of the trajectory 
is such that yl ~ 1, the formation length is very short {Lf ~ A), and the 
radiated field is practically zero. As a result, Eq. (9) can always be applied, 
i.e. the paraxial approximation can always be applied, whenever 7^ » 1. 

Complementarily, it should also be remarked here that the status of the 
paraxial equation Eq. (9) in Synchrotron Radiation theory is different from 
that of the paraxial equation in Physical Optics. In the latter case, the paraxial 
approximation is satisfied only by small observation angles. For example, 
one may think of a setup where a thermal source is studied by an observer 
positioned at a long distance from the source and behind a limiting aperture. 
Only if a small-angle acceptance is considered the paraxial approximation 
can be applied. On the contrary, due to the ultra-relativistic nature of the 
emitting electrons, contributions to the SR field from parts of the trajectory 
with formation length » A (the only non-negligible) are highly colli- 
mated. As a result, the paraxial equation can be applied at any angle of 
interest, because it practically returns zero field at angles where it should 
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not be applied. 



Finally, since the characteristic scale of variation of E is much larger than A, 
the paraxial approximation is valid up to distances of the observer from the 
electromagnetic sources of order A. 

The Green's function for Eq. (9), namely the solution corresponding to the 
unit point source can explicitly be written in an unbounded region as 



G(r- f , z-z') = ; ■ exp 

Note that when z-z' < the paraxial approximation does not hold, and the 
paraxial wave equation Eq. (9) should be substituted, in the space-frequency 
domain, by the more general Helmholtz equation. However, the radiation 
formation length for z-z' < is very short with respect to the case z-z' > 0, 
i.e. there is effectively no radiation for observer positions z - z' < 0. As a 
result, in this paper we will consider only z-z' > 0. The reason why |z - z'| 
is present in Eq. (10) (while z - z' > always) is that, mathematically, the 
Green's function G is actually related to the operator on the left hand side 
of Eq. (9), and not to the whole equation. For example, when dealing with 
wavef ront propagation, one must consider the homogeneous version of Eq. 
(9), and the same Green's function in Eq. (10) can be used, as we will see, to 
propagate the electric field. In this case, propagation can be performed in 
the backward direction as well, i.e. for z-z' <0. 

Note that Eq. (10) automatically include information about the boundary 
condition for the field. In the present case, since we are dealing with un- 
bounded space, the field vanishes at large distance from the sources. Due 
to this fact, the Green's function in Eq. (10) is a scalar function, while in 
general it admits tensorial values. Using the definition of Green's function, 
and carrying out integration over transverse coordinates we obtain 



r' 



tco- 



2c\z - z'\ 



(10) 



E = ^j dz'{^^?x(zOG(r-r'o(z'),z-z') + [V'^G(r-'-r^z-z')].^^^^^^^ 

—oo 

's{z) z 



X exp ICO I 



V c. 



(11) 



where indicates derivative with respect to r'. Explicit substitution of Eq. 
(10) yields the following result 
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z 

£(z,f) = — 5- dz ; 



Vi.{z') _ r- rojz') 
c z-z' 



X exp 



ICO 



r- r' |2 r ^ 1 

h I rfz 

2c(z - z') j 2y2(z)c 



(12) 



where we the choice of integration limits in dz indicate that the electron 
arrives at position z = at time ta = 0. Eq. (12) is valid at any observation 
position z such that the paraxial approximation is valid, i.e. up to distances 
between the observer and the electromagnetic sources comparable with the 
radiation wavelength. One may recognize two terms in Eq. (12). The first in 
Vj_{z') can be traced back to the current term on the right hand side of Eq. 
(7), while the second, in r- roi^')/ corresponds to the gradient term on the 
right hand side of Eq. (7). 

Eq. (12) is used as starting point for numerical codes like SRW. The only 
approximation used is the paraxial approximation. This rules out the pos- 
sibility of using SRW to study the region of applicability of the paraxial 
approximation. However, once the paraxial approximation is granted for 
valid, SRW, or Eq. (12), can be used to investigate the applicability of ER 
theory, which is built within the constraints of the paraxial approximation. 
Note that the evaluation of the field begins with the knowledge of the trajec- 
tory followed by the electron, which is completely generic. In other words, 
one needs to know the electromagnetic sources to evaluate the field at any 
position z down the beamline. 

Alternatively, the knowledge of the far-zone field distribution, i.e. a limit 
of Eq. (12), allows one to specify an algorithm to reconstruct the field in 
the near zone up to distances of the observer from the sources much larger 
than A. An important characteristic of this algorithm is that it works within 
the region of applicability of the paraxial approximation, » 1 only. Such 
algorithm was developed in |15|l to deal with SR problems in full generality, 
and will be used in this paper to develop our ER theory. It follows three 
major steps. 



1. The first step is the characterization of ER emission in the far zone. From 
Eq. (12) follows directly: 
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xexp 
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dz 



ICO 

2c 



2cyl{z) 



(13) 



— * 

where 6 - r/z defines the observation direction, and 6 = \d\. Note that the 
concept of formation length is strictly related to the concept of observation 
angle of interest. In fact, given a certain formation length Lf, and substituting 
it into the phase in z'd^ in Eq. (13), one sees that the integrand starts to be 
highly oscillatory for angles 6 ^ -^A/Lf. Eq. (13) is taken as the starting 
point for our algorithm. 



II. The second step consists in interpreting the far-zone field in Eq. (13) as a 
laser-like beam, generated by one or more virtual sources. These sources are 
not present in reality, because they are located at positions inside the mag- 
netic setup, but they produce the same field as that of the real system. Hence 
the denomination "virtual". A virtual source is similar, in many aspects, to 
the waist of a laser beam and, in our case, exhibits a plane wavefront. It is 
then completely specified, for any given polarization component, by a real- 
valued amplitude distribution of the field, located at a fixed longitudinal 
position. 

Suppose we know the field at a given plane at z, and we want to calculate 
the field at another plane at Zj. In paraxial approximation and in free space, 

the homogeneous version of Eq. (9) holds for the complex envelope E of the 

Fourier transform of the electric field along a fixed polarization component, 

-» 

that is [Vx^ + {2ico/c)dz]E = 0. One has to solve this equation with a given 
initial condition at z, which defines a Cauchy problem. We obtain 

£(z„ i) = -— \dr E{z, f)G (f- r, z, - z) , (14) 



where the integral is performed over the transverse plane and the Green's 
function G in unbounded space is given in Eq. (10). Similarly as before, 

-i 

it is important to remark that since £ is a slowly-varying function with 
respect to the wavelength, one cannot resolve the evolution of the field on 
a longitudinal scale of order of the wavelength within the accuracy of the 
paraxial approximation. In order to do so, the paraxial equation should be 
replaced by the more general Helmholtz equation. Let us now consider the 
limit z — > CX3, with finite ratio r'/z. In this case, the exponential function in 
Eq. (14) can be expanded giving 
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Letting 6 - r'/z we have 



2c 



(z + Zs) 



E (z, 0^ exp 



— r ■ 

c 



(16) 



where the transverse vector r defines a transverse position on the virtual 
source plane at z = z,. Eq. (16) allows to calculate the field at the virtual 
source once the field in the far zone is known. The specification of the vir- 
tual source amounts to the specification of an initial condition for the electric 
field, that is then propagated at any distance. From this viewpoint, iden- 
tification of the position z = Zs with the virtual source position is possible 
independently of the choice of Zg. In other words, like in laser physics, the 
SR field can be propagated starting from any point Zg. However, there are 
choices that are more convenient than others, exactly like in laser physics the 
waist plane is privileged with respect to others. The most convenient choice 
of Zg is the one that allows maximal simplification of the phase contained in 
the far-zone field E(0) with the quadratic phase factor in 0^ in Eq. (16). In 
practical situations of interest it is possible to choose z, in such a way that 
the field at the virtual source exhibits a plane wavefront, exactly as for the 
waist of a laser beam. Finally, in some cases, it is convenient to consider the 

far zone field £(0) as a superposition of different contributions. In this way, 
more than one virtual source can be identified and treated independently, 
provided that different contributions are finally summed together. 



III. The third, and final step, consists in the propagation of the field from 

-» 

the virtual sources in paraxial approximation. Each source £(Zs, r) generates 
the field 

~, _A ICO ( , -* ~ -* 
E{z,r) = - — dr' £(Zs,r')exp 

as follows directly from Eq. (14). 

Here it should be stressed that the inverse field problem, which relies on 
far-field data only, cannot be solved without application of the paraxial 
approximation. 

If the paraxial approximation were not applicable, we should have solved 
the homogeneous version of Eq. (2). Boundary conditions would have been 



ico \r - r'\ 
2c(z - Zs) 



(17) 
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constituted by the knowledge of the field on a open surface (for example, 
a transverse plane) and additionally, by Rayleigh-Sommerfeld radiation 
condition at infinity, separately for all polarization components of the field. 
However, this would not have been enough to reconstruct the field at any 
position in space. In order to do so, we would have had to specify the 
sources. In fact, the boundary conditions specified above allow one to solve 
the direct transmission problem, but not the inverse one. 

Since, however, the paraxial approximation is applicable, the inverse field 
problem turns out to have unique and stable solution. A conservative esti- 
mate of the accuracy of the paraxial approximation is related to the distance 
d between the point where we want to know the field and the electron trajec- 
tory in the space-frequency domain [jist]. When d > Lf this accuracy can be 
estimated to be of order A/Lf.It quickly decreases as Ly^ » » A remaining, 
at least, of order oi A/d. 

The fact that we can reconstruct the near-zone field starting from the knowl- 
edge of the far-zone field and from the propagation equation looks para- 
doxical. In fact, in the far zone, all information about the velocity field in the 
Lienard-Wiechert expressions for the field is lost. It is interesting to note, for 
example, what is reported in [28] concerning the velocity term: "In the case 
of infrared synchrotron radiation, and THz radiation in particular, this term 
is not small and must be included in all calculations". Here we apparently 
seem to have lost track of every information about the velocity term. As 



shown in [ISD, the paradox is solved by the fact that, although in the far- 
zone limit of Eq. (12) includes information about the Fourier Transform of 
the acceleration term of the Lienard-Wiechert fields only, information about 
the velocity term is included in the field propagation equation through the 
Green's function Eq. (10), which solves Maxwell's equations as the Lienard- 
Wiechert expressions do. 



4 Sharp-edge approximation 



Let us consider the system depicted in Fig. 1(a). In this case of study the 
trajectory and, therefore, the space integration in Eq. (13) can be split in 
three parts: the two bends, which will be indicated with bi and b2, and the 
straight section AB. One may write 



£(z, r) = Efoi(z, r) + Eab{z, r) + Ebiiz, r) , (18) 

with obvious meaning of notation. We will denote the length of the segment 
AB with L. This means that points A and B are located at longitudinal 
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coordinates za = 



-L/2 and Zg = L/2. 



Recalling the geometry in Fig. 1(a) we have yziz) = y for Za < z' < Zg. With 
the help of Eq. (13) we write the contribution from the straight line AB 




e'z z' 1 , 



(19) 



where we assumed that ro(z') - 0, i.e. that the particle has zero offset and 
deflection. From Eq. (19) one obtains: 



~ icoeL 
Eab = exp 
c^z 



icod^z 
2c 



6 sine 



4c 



7^ 



(20) 



Note that Eq. (20) describes a spherical wave with the center in the middle 
of the straight section. Moreover, it explicitly depends on L. 

In the previous Section we defined the formation length as the length needed 
for the phase of the electric field seen by an observer on the z axis to overtake 
the phase of the sources of a radiant. It follows from this definition, and from 
the phase in Eq. (19) that the formation length Lf for the straight section AB 
can be written as L^ ~ min[y^A, L]. Then, either Lf ~ y-^A or Lf ~ L. In 
both cases, with the help of the phase in the integrand in Eq. (19) we can 
formulate on a purely mathematical basis an upper limit to the value of the 
observation angle of interest for the straight line, < A/Lf. Moreover, if 
Lf ^ y^A, the maximal angle of interest is independent of the frequency and 
equal to 1/y, in agreement with what has been said before. 

According to the superposition principle, one should sum the contribution 
due to the straight section to that from the bends. However, as discussed in 
Section 2, when 6^1 one can ignore the presence of the bending magnets 

with good accuracy. Note that a direct confirmation of this fact can be given 

-> 

by analyzing explicitly the field from the half bends, e.g. E|,2. An expression 
-> 

for the quantity Eh2 can be found from first principles, applying Eq. (13) to 
the case of a half bend. 



In the framework of the paraxial approximation we obtain for z > L/2: 



s(z) ^ z + 



(z - L/2f 
6R2 



r(s) = -R 



cos 



(^)] 



(21) 



Cx (or gy) being a unit vector along the x (or y) direction. Substitution in Eq. 
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(13) and use of Eq. (8) giveS 



~ icoe 

Ebiiz, ^) = -q exp [zOs] exp [zO, 



oo 
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(22) 



where 



a;z(02 + 02) 




Analysis of the phase term in z'^ Eq. (22) shows that the integrand starts to 
exhibit oscillatory behavior within distances of order of Lf^ = '\/R^A, that is 
the radiation formation length for the bending magnet at A » R/y^- Simi- 
larly, we have seen from Eq. (19), and has been also shown with qualitative 
arguments in Section 2, the formation length for the straight section amounts 
to Lf = min[L, y^A]. The ratio Lf^/Lf is responsible for the relative weight 
of ER compared to bending magnet field contribution. Note that strictly 
speaking, when (p » 1, Lfi,/Lf is equal to 6 ■ (p (and not to 6). However, 
6 ■ (p ^ 1 always, to insure that A » Ac. As a result, in all generality, it is 
possible to talk about ER if and only if 6 « 1 and 6 • « 1 (or '^R^A/L <^ 1, 
A / Ac ^ 1 in terms of dimensional parameters). 

When 6 > 1, one cannot talk about pure ER. One must account for bend- 
ing magnet contributions as well. Then, expressions presented here for the 
electric field from the straight section can be seen as partial contributions, 
to be added to bending magnet contributions calculated elsewhere. 

Note that in this paper we first introduced a measure of "how sharp" the 
edges are through the parameter 6 and, with this, we specified the region of 
applicability for ER theory. 



Usually, in textbooks, the z axis is chosen in such a way that Ox = x/z = 0, i.e. it 
is not fixed, but depends on the observer position. This can always be done, and 
simplifies calculations. However, since the wavefront is not spherical, this way of 
proceeding can hardly help to obtain the phase of the field distribution on a plane 
perpendicular to a fixed z axis. Calculations in our (fixed) coordinate system is more 
complicated and can be found in e.g. in [29]. After some algebraic manipulation 
and change of variables we obtain Eq. (22). 



23 



It should also be noted that bending magnets at the straight section ends 
act like switchers, i.e. they switch on and off radiation seen by an observer. 
Observers see uniform intensity from a bend along the horizontal direction. 
However, not all parts of the trajectory contribute to the radiation seen by 
a given observer, because radiation contributions from different parts of 
the bend is highly coUimated, hence the switching function. Since we are 
not interested in electromagnetic sources responsible for field contributions 
that are not seen by the observer, we may say that bends switch on and off 
electromagnetic sources as well. 

Finally, it should be remarked that the far-zone asymptotic in Eq. (19) is valid 
at observation positions z » L. This is a necessary and sufficient condition 
for the vector n pointing from source to observer, to be considered constant. 
This result is independent of the formation length. When L < y^A we can 
say that an observer is the far zone if and only if it is located many formation 
lengths away from the origin. This is no more correct when L » y^A. In this 
case the observer can be located at a distance z » y^A, i.e. many formation 
lengths away from the origin of the reference system, but still at z ~ L, i.e. 
in the near zone. As we see here, the formation length is often, but not 
always related to the definition of the far (or near) zone. In general, the far 
(or near) zone is related to the characteristic size of the system, in our case 
L. In its turn Lf < L, which includes, when y^A « L, the situation Lf <^L. 

Since in the following we will only deal with a contribution of the electric 

-> 

field, i.e. that from the straight section Eab, from now on, for simplicity, we 
will omit the subscript AB. 

The radiation energy density as a function of angles and frequencies co, i.e. 
the angular spectral flux, can be written as 



rfW _ cz^ 
dcodQ. ~ An^ 



(24) 



dD. being the differential of the solid angle Q. Substituting Eq. (20) in Eq. 
(24) it follows that [iSJ 



dW eWL^ 



dcodQ, 4n^c^ 



02 sinc^ 



4c 



(25) 



It is now straightforward to introduce the same normalized quantities de- 



fined in Section 2: 6 = yfLjA 6, and (p = L/{y^A). We may write the angular 
spectral flux in normalized units as 
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Fig. 6. Normalized angular spectral flux of the radiation from the setup in Fig. 1(a) 
for different values of (p. 



1 = 



coe^ 



(26) 



so that 



2^1 = sinc^ 



(27) 



where z = z/L. Eq. (27) is plotted in Fig. 6 for several values of (/) as a function 
of the normalized angle 6. The natural angular unit is evidently yfXfL. 

The angular spectral flux, once integrated in angles, is divergent, as one can 
see from z^t oc 1/0. This feature is related with the limit of applicability of the 
sharp-edge approximation. Note that within the framework of the paraxial 
approximation alone, the integrated angular spectral flux calculated with 
Eq. (13) does not have any singularity, whatever the electron trajectory 
is. The paraxial approximation, as discussed above, is related to the large 
parameter y^. However, our ER theory is related to another parameter, 
6 1, which controls the accuracy of the sharp-edge approximation. It is 
within the framework of the sharp-edge approximation that the integrated 
flux is logarithmically divergent. Accounting for the presence of the bend 
would simply cancel this divergence. 

We can now justify findings in the previous Section. From Eq. (27) we see 
that, in the limit for (p <^1, the radiation profile is a universal function, and 
peaks at ~ 2.2. When, instead, (/) > 1, radiation is much better coUimated 
peaking at ~ corresponding to ~ 1/y. 
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Fig. 7. Angular spectral flux as a function of the normalized angle 6 for two different 
edge length parameters 6 = 0.2 and 6 = 0.02. Here the straight section length 
parameter (/) ^ 0.01. Left and right plots are obtained cutting the spectral flux 
profile at 0y = and 0^ = respectively (i.e. electron motion is in xz plane). The 
dotted curves are calculated with the analytical formula Eq. (27). Solid lines are the 
results of numerical calculations with computer code SRW. 

The behavior of the far-field emission described here is well-known in lit- 
erature. Nonetheless, the accuracy of the asymptotic expression for 6 « 1, 
Eq. (27), has never been discussed: the parameter 6 has been introduced 
here for the first time. Numerical calculations were never used before to 
scan the parameter space in 6 and to provide an universal algorithm for es- 
timating the accuracy of the ER theory. We can study such accuracy now by 
comparing asymptotical results with SRW outcomes at different values of 
6. This comparison is illustrated in Fig. 7. It can be seen that edge radiation 
approximation provides good accuracy for b < 0.01. 

For completeness, and within the limiting case for 6 «; 1, it is interesting 
to study the accuracy of the asymptotic expression for (/) « 1 of Eq. (27). 
In this case one does not need comparisons with SRW results, because the 
asymptotic limit of Eq. (27) is simply 



zH^e^ sinc^ 



4 



(28) 



Results are shown in Fig. 8. It can be seen that the asymptotic expression 
Eq. (28) provides good accuracy for (p < 0.1. 
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Fig. 8. Angular spectral flux as a function of the normalized angle 6 for different 
straight-section length parameters calculated after Eq. (27) and comparison with 
the asymptotic limit for 1 in Eq. (28). 

From now on, we will only consider the asymptote for sharp-edges 6 ^ lin 
the far-zone. This is the starting point for further investigations of near-zone 
ER, based on the use of virtual source techniques. 



5 Near field Edge Radiation theory 



5.1 Edge Radiation as a field from a single virtual source 



Eq. (16) and Eq. (20) allow one to characterize the virtual source through 
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(29) 



Eq. (29) is valid in any range of parameters, i.e. for any choice of However, 
the Fourier transform in Eq. (29) is difficult to calculate analytically in full 
generality. Simple analytical results can be found in the asymptotic case for 
<^c 1, i.e. for L/(y^A) <^c 1. In this limit, the right hand side of Eq. (29) can 
be calculated with the help of polar coordinates. An analytic expression for 
the field amplitude at the virtual source can then be found and reads: 
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E(0,f) = -,^?smc[-], (30) 

where = |r| as usual. It is useful to remark, for future use, that similarly 
to the far-field emission Eq. (20), also the field in Eq. (30) explicitly depend 
on L. 



It is interesting to comment on the meaning of the phase in Eq. (30), i.e. on 
the factor -i in front of the right hand side. Such phase is linked with the 
(arbitrary) choice of phase of the harmonic of the charge density at z = 0. In 
particular, such phase was chosen to be zero at z = 0. Propagating Eq. (30) 
to the far-zone, one obtains Eq. (20). In other words, the plane wavefront 
transforms into a spherical wavefront centered at z = 0. Note that there is an 
imaginary unit i in front of Eq. (20), meaning that an extra minus sign, i.e. a 
phase shift of n, results from the propagation of Eq. (30). This extra phase 
shift of n represent the analogous of the Guoy phase shift in laser physics, 
and is in agreement with our interpretation of the virtual source in Eq. (30) 
as the waist of a laser-like beam. Note that, while for azimuthal-symmetric 
beams the Guoy phase shift is known to be 7t/2, this result is not valid in our 
case where the cartesian components of the field depend on the azimuthal 
angle. 

We define the normalized transverse position f = r/ VAL. Moreover, since 
the source is positioned at z = 0, we indicate the normalized spectral flux at 
the virtual source as Ig, defined similarly as 1 in Eq. (26). It follows that 



!s{r) = 16 f^sinc^ [f) . (31) 

The profile in Eq. (31) can be detected (aside for scaling factors) by imaging 
the virtual plane with an ideal lens, and is plotted in Fig. 9. 

Note that Eq. (30) describes a virtual source characterized by a plane wave- 
front. Application of the propagation formula, Eq. (14), to Eq. (30) allows 
one to reconstruct the field both in the near and in the far region. We obtain 
the following result: 
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exp 
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(32) 



where we defined 6 - r/z, independently of the value of z. This definition 
makes sense whenever z 0, and yields usual angular distributions in the 
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Fig. 9. Normalized spectral flux at the virtual source. Is, as a function of f (upper 
plot) and 3D view as a function of x and y. 

far zone, for z L. Eq. (32) solves the field propagation problem for both 
the near and the far field in the limit for <^1. 

Eq. (32) is singular at r = (i.e. = 0) and z = L/2. Within our sharp-edge ap- 
proximation, this singularity is mathematically related to the divergence of 
the integrated spectral flux in the far zone, which has been discussed above. 
If one goes beyond the applicability region of the sharp-edge approxima- 
tion by specifying the nature of edges and calculating the field within the 
framework of the paraxial approximation alone (i.e. with Eq. (13)), one sees 
that the integrated angular spectral flux is not divergent anymore, and that 
the field reconstructed at the point r - 0, z = L/2 using this far-field data 
has no singularity at all. 

Note that in the limit for z » L Eq. (32) transforms to 



~ icoeL 
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(33) 
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Fig. 10. Evolution of z^I for edge radiation in the limit for «; 1. These profiles, 

calculated with Eq. (34), are shown as a function of angles at different observation 
distances z = 0.6, z - 1.0, z - 2.0 and z = 5.0 (solid lines). The dashed line always 
refers to the far-zone asymptote, Eq. (28). 

corresponding to the limit of Eq. (20) for L <c y^A, i.e. ^ <sc 1. 

The normalized angular spectral flux associated with Eq. (32) is given by 



V / 02 



exp 



2(1 +2z). 



exp 



2(-l + 2z)- 



(34) 



This notation is particularly suited to discuss near and far zone regions. Eq. 
(34) reduces to Eq. (28) when z » 1. To sum up, when (p -^1 we have only 
two regions of interest with respect to z: 

• Far zone. In the limit for z » 1 one has the far field case, and Eq. (34) 
simplifies to Eq. (28). 

• Near zone. When z < 1 instead, one has the near field case, and Eq. (34) 
must be used. 



Of course, it should be stressed that in the case z < 1 we still hold the 
assumption that the sharp-edge approximation of is satisfied. It is interesting 
to study the evolution of the normalized angular spectral flux for edge 
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(|i = 10.0 



= 50.0 



Fig. 11. Normalized spectral flux at the virtual source for the setup in Fig. 1(a). These 
profiles are shown for = 0.1, (p = 1.0, (p = 10.0 and (p = 50.0 (solid lines). Solid 
curves are calculated with the help of Eq. (35). The dotted lines show comparison 
with the asymptotic limit for <p <^1, shown in Fig. 9 and calculated using Eq. (31). 

radiation along the longitudinal axis. This gives an idea of how good the 
far field approximation (z » 1) is. A comparison between z^I at different 
observation points is plotted in Fig. 10. 

The case (p studied until now corresponds to a short straight section, in 
the sense that L <?c y^A. When this condition is not satisfied, we find that the 
integral on the right hand side of Eq. (29) is difficult to calculate analytically. 
Thus, the single-source method used here is advantageous in the case c/) «; 1 
only. However, one can use numerical techniques to discuss the case for any 
value of (p. With the help of polar coordinates, the right hand side of Eq. (29) 
can be transformed in a one-dimensional integral, namely 



£(0,r-) = 
yielding 



c2 rj 02 + ] 



1/r 



sm 



4c \ 
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dd, 



(35) 
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(36) 



We calculated the spectral flux associated with the virtual source for values 
(p = 0.1, (p = 1, (p = 10 and (p = 50, the same values chosen for Fig. 6. We 
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Fig. 12. Normalized spectral flux at the virtual source for the setup in Fig. 1(a) for 
(j) = 50 (enlargement of the bottom right graph in Fig. 11). 

plot these distributions in Fig. 11. It is also instructive to make a separate, 
enlarged plot of the case (p = 50, that is in the asymptotic case for » 1. 
This is given in Fig. 12. Fine structures are now evident, and are consistent 
with the presence of fine structures in Fig. 6 for the far zone. 

We managed to specify the field at the virtual source by means of numerical 
techniques, even in the case (/) » 1 (see Fig. 12). Once the field at the virtual 
source is specified for any value of (p, Fourier Optics can be used to propagate 
it. However, we prefer to proceed following another path. There is, in fact, 
an alternative way to obtain the solution to the field propagation problem 
valid for any value of (p, and capable of giving a better physical insight for 
large values of (p. 



5.2 Edge Radiation as a superposition of the field from two virtual sources 
Consider the far field in Eq. (20). This can also be written as 




(37) 



where 
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(38) 



When z » L, the two terms £i and E2 represent two spherical waves re- 
spectively centered at Zgi = -L/2 and Zs2 = L/2, that is at the edges between 
straight section and bend^ Analysis of Eq. (38) shows that both contri- 
butions to the total field are peaked at an angle of order 1/y. While the 

-i 

amplitude of the total field explicitly depends on L, the two expressions Ei 

and E2 exhibit dependence on L through phase factors only This fact will 
have interesting consequences, as we will discuss later. The two spherical 

waves represented by Ei and E2 may be thought as originating from two 
separate virtual sources located at the edges between straight section and 
bends. One may then model the system with the help of two separate vir- 
tual sources, and interpret the field at any distance as the superposition 
of the contributions from two edges. It should be clear that contributions 
from these edges are linked with an integral of the trajectory followed by 
the electron. Thus the word "edge" should be considered as a synonym for 
"virtual source" here. 

Let us specify analytically the two virtual sources at the edges of the setup. 
To this purpose, we take advantage of Eq. (16) with Zs = Zs(i,2)/ separately 

substituting Ei and E2 and using polar coordinates. We find the following 
expressions for the field at the virtual source positions Zgi = -L/2 and 
Zs2 = L/2: 



Lsl,: 



,s2 



L \ 2coe 



_ icoL 




(39) 



where Ki{-) is the modified Bessel function of the first order. Analysis of Eq. 
(39) shows a typical scale related to the source dimensions of order Ay in 
dimensional units, corresponding to 1/ -yj^ in normalized units. Also, the 
fact that the field in the far zone, Eq. (20), exhibits dependence on L only 



^ At first glance this statement looks counterintuitive. In order to find where the 

spherical wave is centered, one needs to know where the phase becomes zero. Now, 

-> 

when z = L/2, the phase in 0^ for E2 in Eq. (38) is not zero, hinting to the fact that 

-> 

the spherical wave E2 is not centered at z = L/2. This last observation, however, is 

misleading. In fact, one should account for the fact that the definition 6 - r/z is 

^ ^ _j 

not natural for Ei and £2- In fact, according to it, 6 is measured from the center of 

-» 

the straight section, while it should be measured from z = -L/2 for £1 and from 

-» 

z - L/2 for £2- Also note that Eq. (38) is only valid in the limit L. 
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through phase factors implies that the field at the virtual sources, Eq. (39), 
exhibits dependence on L only through phase factors (and viceversa). 

Application of the propagation formula Eq. (14) allows to calculate the field 
at any distance z in free-space. Of course, Eq. (39) can also be used as input to 
any Fourier code to calculate the field evolution in the presence of whatever 
optical beamline. However, here we restrict ourselves to the free-space case. 
In order to simplify the presentation of the electric field we take advantage 

of polar coordinates and we use the definition E = EyAL c/e (so that 7, 

-> 

introduced in Eq. (26), is given by 7 = |£p) for the field in normalized units. 
We obtain: 
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(40) 



In the limit for z » 1, using Eq. (40) and recalling ^ df f'Ki{ ^f')}i{df') = 
d/[ yf^{6^ + (p)] we obtain back Eq. (27). Similarly, in the limit for (f) ^ 1, 
and using the fact that Ki( y/^f) ^ l/(f y/^) one recovers Eq. (32). In general, 
the integrals in Eq. (40) cannot be calculated analytically, but they can be 
integrated numerically. 



5. 3 Classification of regions of observation 



Qualitatively, we can deal with two limiting cases of the theory, the first 
for «; 1 and the second for » 1. As for the case of a single virtual 
source, there are no constraints, in principle, on the value of (p. However, as 
we will see, the two-source method gives peculiar advantages in the case 
» 1, while, has we have seen before, the case <c 1 is better treated in 
the framework of a single source. 
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5.3.1 Case cp 



Let us briefly discuss the case (p 1 in the framework of the two-source 
method. In this case, one obviously obtains back Eq. (32). An alternative 
derivation has been shown in Section 5.1. As one can see Eq. (32) is indepen- 
dent of (p. In Fig. 10 we plotted results for the propagation according to Eq. 
(34). Radiation profiles are shown as a function of angles 6 at different ob- 
servation distances z = 0.6, z = 1.0, z = 2.0 and z = 5.0. As discussed before, 
one can recognize two observation zones of interest: the near and the far 
zone. As it can be seen from Eq. (32), the total field is given, both in the near 
and in the far zone, by the interference of the virtual source contributions. 
The virtual sources themselves are located at the straight section edges. Eq. 
(40) shows that the transverse dimension of these virtual sources is given by 
yA in dimensional units. This is the typical scale in r' after which the inte- 
grands in df' in Eq. (40) are suppressed by the function Ki. Thus, the sources 
at the edges of the straight section have a dimension that is independent of 
L. In the center of the setup instead, the virtual source has a dimension of 
order VAL as it can be seen Eq. (31). When (p <^ 1 the source in the center 
of the setup is much smaller than those at the edge. This looks paradoxical. 
The explanation is that the two contributions due to edge sources interfere 
in the center of the setup. In particular, when (p <^ 1 they nearly compen- 
sate, as they have opposite sign. As a result of this interference, the single 
virtual source in the center of the setup (and its far-zone counterpart) has a 
dimension dependent on L (in non-normalized units) while for two virtual 
sources at the edges (and in their far-zone counterpart) the dependence on 
L is limited to phase factors only. Due to the fact that edges contributions 
nearly compensate f or (/) «; 1 one may say that the single-source picture is 
particularly natural in the case ^ <^ 1. 



5.3.2 Case » 1 

Let us now discuss the case (/) » 1. In this situation the two-sources picture 
becomes more natural. We indicate with di^z - z ± L/2 the distances of the 
observer from the edges. From Eq. (19) we know that when » 1 the 
formation length is Lf = y^K, much shorter than the system dimension L. 
As a result, one can recognize four regions of observation of interest. 

In Fig. 13 we plotted, in particular, results for the propagation in case (p) — 
50. In this case, for arbitrary z, integrals in Eq. (40) cannot be calculated 
analytically, but they can be integrated numerically. Radiation profiles are 
shown as a function of angles Q at different observation distances z = 0.52, 
z = 0.6, z = 1.5 and z = 100.0. 

• Two-edge radiation, far zone: d\;i » L (i.e. z » L). Eq. (20) and Eq. (25) 
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(t) = 50.0 



Fig. 13. z^I at (/) = 50. These profiles are shown as a function of angles at different 
observation distances z = 0.52, z = 0.6, z = 1.5 and z = 100.0. 



should be used. When di^2 ^ L we are summing far field contributions 
from the two edge sources. This case is well represented in Fig. 13 for 
z = 100, where interference effects between the two edges contribution 

are clearly visible. 

Two-edge radiation, near zone: di^2 ~ Eq. (40) should be used. When 
di^2 ~ the observer is located far away with respect to the formation 
length of the sources. Both contributions from the sources are important, 
but that from the nearest source begins to become the main one, as di and 
d2 become sensibly different. This case is well represented in Fig. 13 for 
z = 1.5. 

Single-edge radiation, far zone: y-^A <^ d2 <^ L and r <^ L/y. Eq. (41) 
should be used. When y-^A <^ d2 ^ L, the contribution due to the near 
edge becomes more and more important. Such tendency is clearly de- 
picted in Fig. 13 for z = 0.6. Interference tends to disappear as the near 
edge becomes the dominant one. In this case, one finds that the electric 
field in Eq. (40) reduces to 
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where ^ = r/(z - L/2) = rjdz. Note that ^ is used here in place of 6, 
because by definition = r/z, where z is calculated from the center of 
the straight section, whereas the definition of ^ is related to the edge 
position at z = L/2. It should be remarked that Eq. (41) constitutes the 
field contribution from the downstream edge of the straight section, and 
that the contribution from the upstream edge (at z = -L/2) can be found 
from Eq. (41) by performing everywhere in Eq. (41), i.e. also in ^, the 
substitution L/2 — > -L/2, and by changing an overall sign. Eq. (41) 
corresponds to an angular spectral flux [^SJ 

It is important to specify the region of applicability of Eq. (41) in the 
transverse direction. For a single edge in the far zone, the amplitude of the 
field decreases as r~, as can be checked by substituting the definition of ^ 
in Eq. (41), and does not depend on z nor y for angles of observation larger 
than 1/y. Since ^2 ^ L, such dependence holds for the upstream edge at 
r > L/y, and for the downstream edge at r ~ ^2/7 ^ L/y. As a result, 
for r > L/y, the contributions from the two edges are comparable, and a 
single-edge asymptote cannot be used. It follows that Eq. (41) applies for 
r « L/y. 

• Single-edge radiation, near zone: 

a/AR2 ^ d2 < y^A, r « yA. Eq. (43) 
should be used. When VaR^ « ^2 ^ y^A we have the contribution from 
a single edge in the near zone. 

As ^2 becomes smaller and smaller the maximum in the radiation profile 
increases (see Fig. 13). This behavior is to be expected. In fact, on the one 
hand the virtual source exhibits a singular behavior at r = 0, while on the 
other hand the integral in Eq. (14) must reproduce the virtual source for 
z — > Zg. In other words, for z — > Zg, the propagator must behave like a 
Dirac 6-distribution. However, the way such asymptote is realized is not 
trivial. At any finite distance dz from the source, Eq. (14) eliminates the 

singularity of the Bessel Ki function. This means that the maximum value 
-» 

in |Ep increases as ^2 decreases, but it always remains finite. In particular, 
-* 

at = 0, |Ep = 0. However, by conservation of energy the integral of 

|£p = over transverse coordinate must diverge at any finite distance 
from the source, because the field diverges at the source position. 

Note that when r ^ yA and di ^ y^A, the integral pertaining the near 
(downstream) edge (at Zs2 = L/2) in Eq. (40) can be calculated analytically. 
In fact, the Bessel Ki function in the integrand can be expanded for small 
values of the argument when f' « 1 / -yj^. When this is not the case 

^ 1 / phase factor under the integral sign makes the integrand 

exhibiting oscillatory behavior (because z - 1/2 «c l/cp, since d2d <sc y^A). 



37 



Contributions to the integrals are therefore negligible. As a result, in this 
case one can use the expansion Ki{ -yj^f) ~ l/( ■\f^f'). Then, using the fact 
that r dxJiiAx) exp[iBr^] = l/A{l-exp[-zAV(4B)]} (for A and B positive). 



one obtains 
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(43) 



Note that while the modulus of Eq. (43) is independent of (p, its region of 
applicability is related to (p and the asymptotic expression deviates from 
Eq. (40) for smaller value of 6 when cp is larger. In fact, Eq. (43) is valid 
only when z - 1/2 <^ l/(p and z0 « 1/ (i.e. r ^yA and ^2 ^ y^A). 

It is interesting to remark here that Eq. (32), which was derived for 
(p ^1, reduces to Eq. (43) when d2 <^ L and r « VAL. It follows that the 
validity of Eq. (43) has a wider region of applicability than that considered 
here. In fact, it may be applied whenever 's/AB^ <^ d2 <^ min(L, y^A) and 

r «c A/Ainin(L~7^. 

If we propagate Eq. (43) to the far zone we obtain an asymptote which 
is valid only for angles much larger than 1/y (i.e. Eq. (43) is an asymptote 
for high values of spatial frequencies). The modulus of Eq. (43) does 
not depend on y (while in the non-asymptotic case radiation for any 
value of z must depend on y, because the far-field radiation from a single 
edge depends on y too) nor it includes information about distribution in 
the far zone within angles comparable with 1 /y. Thus, the applicability 
of this high-spatial frequency asymptote depends on what practical (or 
theoretical) problem we try to solve. It is useful, for example, if we discuss 
about a sample in the very near zone. However, if we discuss about design 
of beam line with an acceptance angle comparable with 1/y (which is 
equivalent to some spatial-frequency filter) the asymptotic expression in 
Eq. (43) cannot be applied anymore, and one should use exact results from 
the propagation integral, i.e. the near-field expression Eq. (40). 

The above-given classification in zones of interest with asymptotical expres- 
sions for the electric field constitutes an important result of our paper. In fact, 
expressions for the electric field without explicit specification of their region 
of applicability are incomplete, and have no practical nor theoretical utility. 
From this viewpoint, it is interesting to compare our results with literature. 
We will limit our discussion to a comparison with recent review | |28|l , which 
summarizes up-to-date understanding of ER within the SR community. 



One result in [281] (Eq. (26)) corresponds to the square modulus of our Eq. 
(43), the single edge near-zone case. The region of applicability specified 
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in [|28|l for such resul10 is A ^ ^2 ^ y^A and L — > oo. A first problem in 
applying this prescription is intrinsic in the condition L — > oo, as there is 
no comparison of L to any other characteristic length. Secondly, this result 
is independent of y and L. As a result, it cannot be valid for arbitrary 
transverse distance r. In fact, since the far-zone field depends on both y and 
L, should we propagate Eq. (43) in the far zone, we could never obtain an 
outcome dependent on y and L. It follows that, as we discussed above, Eq. 
(43) is valid for arbitrary L (under the sharp-edge limit ^/AF? «: L), but its 
region of applicability depends on L or y: '\/AR^ « ^2 ^ min(L,y^A) and 
r « -yjA min(L, y^A). Note that the requirement r « Ay for the applicability 
of Eq. (43) is present in the or igin al paper [8], but has been omitted in 
some later publications, e.g. |@,[lO, isj]. Because of this the dependence of 



near-zone single-edge radiation upon y was unclear, and the asymptotic 
expression Eq. (43) started to be considered by other authors, as in [28], 
without proper requirements on r. 

A second result in [28] for the case of a finite straight section length L (Eq. 
(27)) corresponds to our Eq. (34). This fact can be proved with the help of 
straightforward mathematical steps. The region of applicability of Eq. (27) 
in |28] is specified by the words "under conditions of validity of equation 
(26)", i.e. A « ^2 ^ y^A. In contrast to this, we have seen that Eq. (34) is 
valid for arbitrary ^2 ^ '^AR^, but includes limitations on L in the form: 
i/A^ « L « y^A. 

The following result in 11281], Eq. (28), is the far-zone single-edge result, i.e. 
our Eq. (42). Eq. (28) in |28l] is presented without region of applicability 
while, as we have seen, it is valid for y^A <^ d2 <^ L and r <^c L/y. Note 
that in this case the denomination "far zone" is not related with the usual 
understanding ^2 — > as is the case for the usual far-zone expression for 
two-edges, but it includes a limit on d2, i.e. ^2 ^ ^■ 



The final result in |28|1 is Eq. (29), that is our far-zone two-edge case, Eq. (25). 
Also Eq. (29) is presented without region of applicability while, as we have 
seen, ^2 ^ I- 

It should be appreciated how our analysis of ER through the parameter 6 
allowed us to define "how sharp" the edges are, and to specify the region 
of applicability for ER theory. 

As a final remark, note that in literature the single edge far-zone case is 
usually presented as the simplest and fundamental case, while in our view 
this is the most complicated and misleading case to discuss. For the sake 
of exemplification, consider the usual assumption made in this case, i.e. 



Converted to our notation. 
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Fig. 14. Transition undulator radiation geometry. 

L — > oo. In order to discuss the far-zone expression, one needs z — > oo, and 
comparing z to L becomes impossible. In other words, L drops out of from 
the problem parameters and, as in [28], it never appears in the condition for 
the region of applicability anymore. In contrast to this, our simplest model 
is the two-edge far-zone model, whose region of applicability is: z » L (or 
^2 ^ It is independent of y and it is much easier (although an extra- 
limitation on angles « i/A/R should be included, due to the sharp-edge 
approximation). After introduction of this model, the following natural step 
was to generalize it to the near zone, introducing more complicated regions 
of applicability discussed before. 



6 Transition undulator radiation 



In this Section we apply the method of virtual sources to the more compli- 
cated case of an undulator setup. 

Instead of the setup in Fig. 1(a), we now consider the system depicted in Fig. 
14 and we consider a single particle moving along the system. The electron 
enters the setup via a bending magnet, passes through a straight section 
(segment AB), an undulator (segment BC), and another straight section 
(segment CD). Finally, it leaves the setup via another bend. Radiation is 
collected at a distance z from the center of the reference system, located in 
the middle of the undulator. The passage of the electron through the setup 
results in coUimated emission of radiation in the range A » A,, and A » Ac, 
where A^ is the resonance wavelength of the fundamental harmonic of the 
undulator, i.e. the extra characteristic length introduced in the setup. This 
kind of radiation is known in literature as Transition Undulator Radiation 
(TUR) iBHBHHl- We will retain this name although, as we will see. 



what we are really discussing about is edge radiation from an undulator 
setup. 
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In our case of study the trajectory and, therefore, the space integration in 
Eq. (13) can be split in five parts: the two bends, which will be indicated 
with bi and bz, the two straight sections AB and CD and the undulator BC. 
One may write 

^ ^ ^ ^ ^ ^ 

E{z, r) = Ebi{z, r) + Eab{z, r) + £bc(Z/ r) + £cd(Z/ r) + Ei,2{z, r) , 

(44) 

with obvious meaning of notation. 

We will denote the length of the segment AD with Ljot/ while we will indicate 
the length of the straight section AB with Li, the length of the straight 
section CD with L2 and the length of the undulator with L^u,. It follows 
that Ltot = Li + Lu, + L2. This means that point A is located at longitudinal 
coordinate Za = -Li - La,/2, while B, C and D are located respectively at 
Zb = -Lu,/2, zc - Lu,/2 and Zo - L^,/! + L2. 

We will describe the field from our TUR setup as a superposition of three 
laser-like beam from straight sections and undulator. As before, with the 
help of Eq. (13) we will first derive an expression for the field in the far zone. 
Then we will calculate the field distribution at the virtual source with the 
help of Eq. (16). Finally, Eq. (17) will allow us to find an expression for the 
field both in the near and in the far zone. 

6.1 Far field from the undulator setup 

Let us describe the far field from the undulator setup in Fig. 14 by separately 
characterizing different field contributions and finally adding them together. 

6.2.2 Field contribution calculated along the undulator 

-> 

We first consider the contribution Ebc from the undulator. Assuming a pla- 
nar undulator with N^, periods we write the following expression for the 
transverse velocity of an electron: 

cK 

v_i{z) = sin {ku,z)x . (45) 



Here K = (Ai^eHiy) l{2nm^c^) is the undulator parameter, being the electron 
mass and Hy, being the maximum of the magnetic field produced by the 
undulator on the z axis. Moreover, k^u, = 2n/Ajv, where Ajj, is the undulator 



41 



period, so that the undulator length is Lu, = NjyAa,. The transverse position 
of the electron is 

n)(z) = —r- cos {k^z)x . (46) 

ykw 



We can now substitute Eq. (46) and Eq. (45) in Eq. (13). Such substitution 
leads to an expression valid in the far zone. We obtain 



~ . ^ . icoe r , , r-^ 1 
EBc{z,r,co) = -^ I dz exp [zOfic] 



-sin(fcj(,z') + 



f + dyg\ . (47) 



Here 



where the average longitudinal Lorentz factor fz in Eq. (48) turns out to 
be 7z = y/ yjl + K^ll, and is always smaller than y because the average 
longitudinal velocity of the electron inside the undulator is smaller than 
that along the straight sections. 

In this paper we will be interested up to frequencies much lower than the 
resonance frequency, i.e. A » A^-, with A^ = ll{l'Yycw)- 

We can show that this condition is analogous, for TUR radiation, to condition 
• 6 «; 1 f or the simple edge radiation setup in Fig. 1(a). In order to do so, 
we first need to discuss the formation length associated with the undulator, 
i.e. with Eq. (47). The definition of formation length was introduced before 
as the value of Z2 - Zi for which the right hand side of Eq. (8) is of order 
unity. However, the physical meaning of formation length is related with 
the integration range in z' such that the integrand in Eq. (47) exhibits an 
oscillatory behavior. In our case, not only the phase in Eq. (48), but also 
the sin(-) term in Eq. (47) have an oscillatory character, and must be taken 
into account when calculating the formation length. As a result, in the long- 
wavelength asymptotic (A » A^), we deal with a situation where the sin(-) 

and the terms in Eq. (47) have different formation lengths. From the sin( ) 
term in Eq. (47) and from Eq. (48) follows a formation length L/ = Ai„, and a 
characteristic angle 6 ~ yfXfK^. However, the TUR contribution is given by 
the terms in 6^ and Oy in Eq. (47). For these terms, from Eq. (48) follows a 
formation length ~ min(y^ A, Na;Ai„), in the limit for A » A^. In this limit, 
the TUR contribution is coUimated to angles 6^ <c A/Aa,. Only within these 
conditions one can properly talk about TUR. Note that L/ » A„; always. 
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because ylA ~ Ai„A/A^ » Aa, and we assume » 1- The expression for the 
formation length Lf above is analogous to that for edge radiation, which is 
given by mm{y^A, L). The analogous of the (p parameter is now given by (pj^ = 
^w^wl {yl.^)f while the analogous of the 6 parameter is = A^,/ (Aa,Na,) = 
1/Nk, «; 1, as Na, » 1. It follows that b„ = Aw/ifl^) ~ ^r/A «: 1. 

For (/)a, • 6a, «: 1 and 6zv ^ 1 the contribution due to the term in sin (/Cj^z') 
in Eq. (47) can always be neglected when compared with the maximal field 
magnitude of the terms in dx,y. Similarly, in Eq. (48), phase terms in cos{kwZ') 
and sin(2fca,z') can also be neglected. As a result, Eq. (47) can be simplified 
as 



£bc(Z/ = J dz'exp U^bc] {dxX + dyfj 



(49) 



where 



Obc = - 



(50) 



6.1.2 Fie/rf contribution calculated along the straight sections 

With the help of Eq. (13) we write the contribution from the straight line AB 
as 



ZB 

= ^ J" dz' exp [zO^b] (0;,x + dyy) 



(51) 



za 



where O^b in Eq. (51) is given by 



Oab = - 



02 z' / 1 ^ 

2 2 



2 1 '^w _|_ £-10 

4^1 4^ 



(52) 



The contribution from the straight section CD is similar to that from the 
straight section AB and reads 



ZD 

= ^ J" dz' exp [zOcd] [ej + 0yy) 



(53) 



Zc 
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where Ocd ut. Eq. (53) is given by 



OcD = - 



+ — ^ + + 



(54) 



In general, the phases Ocd and Oab start exhibiting oscillatory behavior 
when z'/(2y^A) ~ 1, which gives a maximal integration range in the longi- 
tudinal direction. Similarly as before, in general one has that the formation 
lengths Lfsi and Lfs2 for the straight sections AB and CD can be written as 

Us(\,2) ~min[A7^L(i,2)]. 

6.1.3 Total field and energy spectrum of radiation 

The contributions for segment AB and segment CD are given by Eq. (51) 
and Eq. (53). One obtains 



— icoeLx 
^AB = exp 



xexp 




(55) 



Similarly, 



~ ia)eL2 
EcD = exp 



c^z 



xexp 



icod^z 
2c 



6 sine 
exp 



4c \y2 



y2 



4c 



(56) 



Finally, the contribution for the segment EC is obtained from Eq. (49). Cal- 
culations yield: 



~ icoeL^, 
Ebc = 



icoO^z 
2c 



6 sine 




(57) 



The total field produced by the setup is obtained by summing up Eq. (55), Eq. 
(56) and Eq. (57). By this, we are neglecting bending magnet contributions. 
A sufficient condition (in addition to the already accepted ones, Nn, » 1 
and Kl Kr » 1) is A » ^/^z- ^^^t, we may neglect bending magnet 
contributions for 61,2 ^ 1, but in this setup Li 2 may be set to zero, in which 
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case we should also impose that the formation length of the bend be much 
shorter than flA, which reduces to A » R/fl- 

As before, the observation angle is measured starting from the center of the 

undulator, located at z = 0, i.e. 6 = r/z. The angular spectral flux can be 
written substituting the resultant total field in Eq. (24). We obtain 



dW 
dcodCl 

+ exp 



. ojLi 
ia)L-a 

+exp 



(1 



+ 



ia)L,n 



exp 

.coU 



.coL 



2/32 



4cy2 \^ 2 ^ 



2/32 



exp 
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4cy2 
, C0L2 



4c> 

1 + 



+ exp 



.COL; 



2a2] 



(i + y'e') 



+ 



icoU,, 



4cy 

2a2 



2il + y + y'0 



'2cy2 V / 4cy2 

that is equivalent to the analogous expression in f^. 



(58) 



Note that Li, L2 and can assume different values, y and y^ are also dif- 
ferent. It may therefore seem convenient to introduce different normalized 
quantities, referring to the undulator and the straight lines. However, in 
the end we are interested in summing up all contributions from different 
sources, so that it is important to keep a common definition of vertical 
displacement (or observation angle). Therefore we prescribe the same nor- 
malization for all quantities: 



a _ /^tot^ , _ hot 



and f = 



(59) 



Then, we introduce parameters Li = Li/Ltot, L2 = L2/Ltot, = L^y/Ltot, 
(/)i,2 = Li^2/(y^A) = Li^2(pt and (pzv = Uv/ifl^), as seen above. Here it should 
be clear that (pt has been introduced only for notational convenience, while 
real parameters related to the physics of the problem are (^1,2- Finally, we 
define Zs = Zs/i-tot- From Eq. (58) follows 



z I- 



402 



{(pt + e^r 
(pf + e^ 



exp 



-■^{(pt + e')-j{(p,. + Lj') 



exp 



-j((P,,+Lj^) 



+ 



(p^/U + 0^ 



exp 



^((p,u + K0^) +exp -j((p,„ + Lu,6^) 



I 

L 4 



45 



exp 



+ exp 



(60) 



where I = E and E = E VALtot c/e. Note that outside the undulator the 

longitudinal velocity is nearer to c than inside {fl < y^). It follows that 
the contribution of the undulator is suppressed compared with that of the 
straight sections, and in the case of comparable lengths and » 1, the 
straight section contribution becomes dominant. 

In the special case for Li = L2 = Lj„ = L/3, (pi^2 = (pt/^, and Eq. (60) simplifies 
to 



z I- 



4:6^ 



-2i- 



4z cos 



-(2d' 



+ 



+ x 



3(/),„/L + 02 



sm 



sin 









(61) 



Eq. (61) can be readily evaluated. As an example, we can calculate the 
intensity distribution of TUR emitted by the SASE 1 European XFEL setup 
at a wavelength A = 400 nm. We assume that the XFEL operates at 17.5 
GeV. Setup parameters are Li = Lj„ = L2 = 200 m, = 400 m, K = 3.3 and 
= 3.56 cm iQ. In this case b ~ 10"^ and ~ 1/N,„ ~ 10"^ while (/)f ^ 8.0 
and (pio ^ 52. Results are plotted in Fig. 15. In that figure, we also propose a 
comparison with outcomes from SRW at z = 6000 m (vertical and horizontal 
cuts). 

In the far zone, well-accepted expressions for the TUR emission are reported 
in literature |16, l^r iM 20], that are equivalent to the following equation for 
the radiation energy density as a function of angle and frequency: 



rfW _ e' 
dcodO. n^c 



(1 + K2/2 + y202)(l + ^202) 



sm 



2y2A 



2n2] 



(62) 



We will show that Eq. (62) cannot be applied for TUR calculations. In our 
understanding, there cannot be any range of parameters in the setup in Fig. 
14 where Eq. (62) is valid. 

In order to prove this it is sufficient to compare Eq. (62) with Eq. (58). Eq. 
(62) does not depend on the straight section lengths Li or L2, and can be 
applied when Li, L2 — > 0. 



However, in that limit, Eq. (58) reduces to: 
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Fig. 15. Cross-check of Eq. (61) with the help of SRW. Here Li = = L2 = 200 
m, E = 17.5 GeV, A = 400 nm, R = 400 m, and K = 3.3. Here A^a = 3.56 cm. The 
observer is located at z = 6000 m. Circles represent horizontal and vertical cuts of 
the intensity profiles calculated numerically with SRW. The solid curve is calculated 
with Eq. (61). The dashed curve is obtained with the near-zone expressions Eqs. 
(69)-(72). 



dW 
dcodO. 



-.2 



1 + + y202 



sm 



(63) 



that is obviously different from Eq. (62). Also note that in the limit for K — 
Eq. (62) tends to zero, whereas Eq. (63) gives back Eq. (25) as it must be. 



6.2 Virtual source characterization and field propagation 



There is a general need, in the EEL community, to extend the current theory 
of TUR to cover the near zone. Eor instance, a possible use of coherent TUR to 
produce visible light synchronized w^ith X-rays from an X-ray free-electron 



laser is discussed in [241 ] . As we have seen, TUR can be discussed as a more 
complicated edge-radiation setup. Within the sharp-edge approximation we 
have contributions from three parts, two straight lines and the undulator. 
The undulator contribution is similar to a straight line contribution, the only 
difference being a different average longitudinal velocity of the electron. 
Then, the far-zone region can be identified by distances z » L^f. 

Expressions in Eq. (55), Eq. (56) and Eq. (57) can be interpreted as far field 
radiation from separate virtual sources. For each far field contribution we 
use a picture with two virtual sources, located at the ends of the straight 
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sections and of the undulator. This makes a total of six sources. However, 
since the virtual source at the downstream [upstream] edge of the first 
[second] straight section has the same longitudinal position of the virtual 
source at the upstream [downstream] edge of the undulator, i.e. z = -La,/2 
[z = La,/2], we combine them together, summing them up by superposition 
principle. As a result, we are left with only four sources, located at 

T Lip a\ 

Zsi = — 2" ~ ~ — 1 ' ~2 ' ^ = — + L2 . (64) 

We obtain an explicit expression for these sources with the help of Eq. (55), 
Eq. (56) and Eq. (57), proceeding analogously as in Section 5.2: 

r \cy I 
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Icoe 



exp 



— exp 



Acyl 



4:cy: 




(66) 
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(67) 



£s4^Y+L2,rj = -— exp 



icoL2 
2cf 



exp 



icoU, 



^cy 




(68) 



In order to calculate the field at any distance z we proceed in analogy with 
Eq. (40), applying the propagation formula Eq. (14). As before, the above 
given equations for the sources can also be used as input to any Fourier 
code to calculate the field evolution in the presence of whatever optical 
beamline. However, here we restrict ourselves to the free-space case. In 
order to simplify the presentation of the electric field we take advantage 

of polar coordinates and we use the definition E = E VALtot c/e (so that /, 
introduced in Eq. (26), is given by 7 = |£p) for the field in normalized units. 
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Note that here z = z/Ltot- We obtain four field contributions, one for each 
source: 



z£i 



111 yf^texp[-iLi(t)t/2] exp [-i(j)J4\ 



e 



z + Lu,/2 + Li 
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2 (z + + Li) 
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(69) 
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(70) 
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(71) 
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2^4 (z, 9 



"^Iz^J^t exp [iL2(ptl^ exp /4] 
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\z — — L21 



X exp 
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i'2 
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exp 



2i>2 



2(z-U/2-L2) 



(72) 



Eqs. (69) to (72) can be used to calculate the field, and hence the intensity, 
at any position of interest in the far and in the near zone. Obviously, in 
the far zone, for z » 1, the square modulus of their sum reduces to Eq. 
(61). As before. It is interesting to cross-check Eqs. (69) to (72) with the 
computer code SRW. We used the same numerical parameters as before: 
Li = = Li ^ 200 m, £ = 17.5 GeV, A = 400 nm and R = 400 m, with 
an undulator parameter K = 3.3. Additionally, we chose different values 
z = 360m, z = 600m, z = 1200m (see Figs. 16 ^ 18), and z = 6000m (see 
Fig. 15), corresponding to z = 0.6, z = 1.0, z = 1.5, z = 2.0 and z = 10 
(which is the far-zone case treated before). Here the bending plane is the 
horizontal plane. As the observation point becomes nearer to the edge of 
the magnet, the influence of the bending magnet becomes more and more 
important, an effect which is evident in the figures from the horizontal cuts 
of SRW two-dimensional intensity profiles. From Figs. 15-rl8 we can see 
that the analytical result for the vertical cut is valid with good accuracy up 
to z = 0.6. 



7 Edge radiation in a waveguide 



In the previous Sections we considered edge radiation propagating in un- 
bounded space. In this Section we extend previous considerations to account 
for the presence of metallic surroundings (like e.g. the electron vacuum 
chamber), which effectively act like a waveguide. 

The need to account for guiding structures in edge radiation setups arises 
when one wants to operate in the sub-millimeter wavelength range. Emis- 
sion of edge radiation in the presence of metallic boundaries has been a 
much less-treated subject in literature, compared to the unbounded space 
case. To the best of our knowledge, there is only one article reporting on 
edge radiation from electrons in a planar overmoded metallic waveguide 
| [l4|] . The analysis in [14] is based on Lienard-Wiechert fields, modified by 
the presence of finite metallic boundary. To account for perfectly conducting 



plates, the author of ||14l] uses a generalization of the method of images, well- 
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Fig. 16. Cross-check of analytical results with SRW. Here Li = L^, = L2 = 200 m, 
E = 17.5 GeV, A = 400 nm, R = 400 m, and K = 3.3. The observer is located at 
z = 360 m. Here A^, = 3.56 cm. Horizontal and vertical cuts of the intensity profiles 
are compared with results obtained with the help of Eqs. (69) to (72). 




Fig. 17. Cross-check of analytical results with SRW. Here Li = La, = L2 = 200 m, 
E = 17.5 GeV, A = 400 nm, R = 400 m, and K = 3.3. Here Ai„ = 3.56 cm. The observer 
is located at z = 600 m. Horizontal and vertical cuts of the intensity profiles are 
compared with results obtained with the help of Eqs. (69) to (72). 
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Fig. 18. Cross-check of analytical results with SRW. Here Li = L^, = L2 = 200 m, 
E = 17.5 GeV, A = 400 nm, R = 400 m, and K = 3.3. Here A-,u = 3.56 cm. The observer 
is located at z = 1200 m. Horizontal and vertical cuts of the intensity profiles are 
compared with results obtained with the help of Eqs. (69) to (72). 

known from electrostatics [31]. Instead of this approach, here we prefer to 
use a mode-expansion technique. We characterize radiation from particles 
in guiding structures with the help of the proper tensor Green's function, au- 
tomatically accounting for boundary conditions. The Green's function itself 
and the field are then conveniently expressed in terms of the natural modes 
of the guiding structure. This approach is shown to provide a convenient 
and methodical way to deal with the vectorial character of our problem. We 
apply our method to the case of a homogeneous waveguide with circular 
cross-section. This configuration approximates the vacuum chamber to be 
used at XFEL and Energy Recovery Linac (ERL) facilities. Our approach, 
however, is very general and may be applied to other geometries, e.g. to 
rectangular waveguides, which approximate vacuum chambers used at SR 
facilities. 



7.1 Definition of setup and position of the problem 



To fix ideas we focus our attention on the setup in Fig. 19. Electrons travel 
through the usual edge radiation setup, similarly as in Fig. 1(a). The differ- 
ence is that now we account for the presence of a waveguide with typical 
cross-section dimension p along the straight section. Since electrons are bent 
away, one may assume that the vacuum chamber ends at the downstream 
bend position z = L/2. Our goal in this Section is to calculate the field dis- 
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Fig. 19. Setup for edge radiation emission includes the presence of a waveguide 
with circular cross-section of radius p along the straight section. 

tribution at that position, which should be subsequently propagated up to 
the experimental station. 

Under the assumption 6 «; 1, which is always considered valid, one can then 
extend the ER theory from unbounded space to the case when a waveguide 
is present. In unbounded space, the assumption 6 «: 1 allowed us to apply 
the ER theory with given constraints on the angular region of observation: in 
fact, at large angles the SR contributions from bends is dominant compared 
to the ER contribution to the field. In other words, as we have seen, ER is a 
more collimated form of radiation, compared to bending magnet radiation. 
The presence of a waveguide seems at first glance to radically modify this 
viewpoint, because both SR from bends and ER is collected and guided into 
the vacuum chamber. However, we will see that this situation still carries a 
strong analogy with the unbounded space case: while in unbounded space 
ER theory can be applied within a given angular region dependent on 6, in 
the waveguide case ER theory can be applied within a given wave-number 
range dependent on 6. As said above, after the field distribution is calculated 
at z = L/2, it has to be propagated up to the experimental station. Radiation 
will be collected within a finite acceptance angle. This finite acceptance can 
be seen as a low-pass filter for spatial frequencies. If the acceptance angle is 
small enough (and this depends on the value of 6) only spatial frequencies 
where ER is dominant will be accepted. Thus, one can conclude that both 
in the case of unbounded space and vacuum chamber, ER theory can be 
applied with given constraints on the acceptance angle of the setup. 

Despite the fact that equations for the edge radiation source are significantly 
complicated by the presence of a waveguide, we find that these complica- 
tions result in the appearance of a single dimensionless extra-parameter. 
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namely the waveguide diffraction parameter Q = p^/{AL). When ^ 1 and 
6 « 1, Q is the only parameter of the problem and has a simple interpreta- 
tion as the ratio between the waveguide cross-section area and the square of 
the radiation diffraction size of the ER virtual source located in the middle 
of the straight section in unbounded space. 

The general theoretical framework of our extended ER theory relies on the 
solution of Maxwell's equations with modified boundary conditions, and 



follows closely the approach proposed in |32l] . We wish to describe, in the 
space-frequency domain, radiation from an ultrarelativistic electron moving 
on a given trajectory inside a metallic vacuum chamber of arbitrary cross- 
section, homogeneous along the z-axis. An example is depicted in Fig. 20, 
where n is a vector field defined on the boundary surface S, such that |n| = 1. 
At any point, n is orthogonal to S and points outwards and at this stage the 
waveguide cross-section is still considered arbitrary. 

Since the electron is ultrarelativistic, the Lorentz factor y obeys <^c 1, 
and paraxial Maxwell's equations can be used in unbounded space. The 
presence of a waveguide introduces the extra-limitation A p. Paraxial 
approximation can be used under these two constraints. Then, the field 
equations read 



^V/E + 2M<9^E) = ^exp 
(hxeX =0 



Is 



(73) 

(vx-e)| =0, 



where the second and the third equation define boundary conditions at the 
waveguide surface S. 



As has been demonstrated in ||33l,l34j,l32l] the field E can be written in terms 



of cartesian components as (compare with Eq. (11)) 
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where derivatives d'^ are taken with respect to r', and -oo has to be taken 
as the point where the source harmonics begin to exist (in our case, at 
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Fig. 20. Geometry of the problem. The electron is moving inside the vacuum cham- 
ber. 

z = -L/2). Eq. (74) is the solution of the inhomogeneous jfield equation, 
describing the jfield from a single electron. As before, ro(z') and Ux(z') fix 
transverse position and velocity of the electron as a function of z', and yz{z) = 
[1 - Vz{z)^/c^]~^^^, Vz indicating the longitudinal velocity of the electron. 
Finally, are the cartesian components of a suitable tensorial Green's 
function for a homogeneous waveguide (along the z-axis): 



icX 
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exp 
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(75) 



where notation ® indicates the tensor product. Here, Cx and Cy are unit vectors 
along cartesian axes x and y, while and A J^'™ are eigenf unctions and 

eigenvalues referring to the scalar problem: 



Vit/;J='™(r) + A™ Vjif) = 
n-(y^ipf)^ =0 

Is 



(76) 



TE and TM respectively standing for Transverse Electric and Transverse 
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Magnetic, with normalization condition 



^ = 1 . (77) 



Note that here we derived the Green's function Eq. (75) directly in paraxial 
approximation starting from Eq. (73) for the transverse components of the 
field. 

Simplifications to Eq. (74) apply remembering that in our case v±{z') - 
and ro(z') - 0. Setting z - L/2 one obtains: 
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(78) 



In a waveguide that is homogenous along the z-axis, only TM modes turn 
out to be driven by the uniform motion of the space-charge distribution. 
As a result, we can focus on TM modes only. This may be directly seen by 
calculating d'^G^. A compact mathematical way proving this fact without 
cumbersome calculations is to note that TE modes appear in the expression 
for the tensorial Green's function, Eq. (75), through the combination V x 
{ip^^ez) = Cx dy i/^J^ - gy dxij^Y' Then, TE modes contribute in d'^ for 

V • [V X (i/^J^^)], which is the divergence of a curl and yields back zero. From 
a physical viewpoint this is a sound result. In fact, radiation is related with 

energy change of the particle, which happens through the scalar product 

-» 

E • V. Since Vj_{z') — 0, TE modes cannot lead to any energy change of the 
electron. As a result, they are not excited. Note that v^{z') = means that 
ro(z') = constant. Since the electron beam trajectory should be centered, one 
recovers ro(z') = 0, leading to extra-simplifications in Eq. (78). 



7.2 Theory of Edge Radiation in a circular waveguide 



We now want to apply Eq. (78) to study the case of the setup in Fig. 19 in the 
presence of a circular waveguide with radius p. Therefore, we fix z = L/2 
in Eq. (78). This allows us to calculate a field distribution at the end of the 
straight section, which acts as a source to be used for further propagation 
in unbounded space or through the beamline optics. 

In order to explicitly compute the tensor Green's function in Eq. (75) we 
need to solve the problem in (76) with this particular choice of geometry. 
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Because of the particular symmetry in the case of a circular waveguide, it 
is convenient to use transverse polar coordinates so that r is identified by 
radial distance r = -^x^ + and angle (p = arctan(i//x). 

Solutions of the equation set in (76) are 



rZ = K^Jm {vnr^r ) sin(m(^) , = A]^U (v„,,r ) cos(m(^) , (79) 

where v,„jt are defined as the roots of /,„ (v^it) = 0, and the normalization 
coefficients for TM modes, A™, are given by 
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(80) 



This leads to the following final result for the tensor Green's function com- 
ponents in cartesian coordinates: 
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where we already dropped TE modes. Eq. (81) is valid in the case of a circular 
waveguide. These facts are demonstrated in |32l] . As a result we have 
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along ix and 
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along gy. A major simplification arises since ro(z') = 0. Because of this, the 
only non-zero contribution to Eq. (82) and Eq. (83) is for m = 0, otherwise 
Jm{Vmkr'lp) yields zero. 

Substitution in Eq. (78) yields, therefore, the following result for the field at 
2 = L/2: 
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Here we defined Q = VQ^^c/(2(i'p^) and = sin((p)^^ + cos{(p)ey. Eq. (84) and 
Eq. (85) solve the problem of characterization of ER in the setup of Fig. 19 
for 6 «; 1. 



It is interesting to show that the treatment of edge radiation in terms of 
virtual sources has a straightforward generalization for the case of a waveg- 
uide. 

It is convenient to introduce the subject by first noting that Eq. (73), the non- 
homogenous Maxwell's equation with boundary conditions at the waveg- 
uide walls, admits a self-reproducing solution. Because of the circular cross 
section of the waveguide, there is no preferential direction on the transverse 
plane. As a result, the transverse electric field must be radially polarized, 

and axially symmetric. We make the ansatz (i) £ - Er{r)Ez{z)er, and (ii) 
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Ez{z) = exp[zz/(2y^A)]. Note that (i) is already equivalent to postulating 
a self-reproducing field, meaning that the transverse shape of the field is 
independent of the longitudinal coordinate z down the waveguide. Substi- 
tuting (i) and (ii) in Eq. (73) and shifting to polar coordinates we obtain the 
following problem for Er{f): 
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where we introduced f = r/ (yA) and p = p/ (yA) ad hoc to simplify notations 
in Eq. (86). Solution of Eq. (86) yields the final result 
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Eq. (87) is well-known (see, e.g. [SSI]), and explicitly describes the self- 
reproducing nature of the field, because the transverse structure of the 
radiation field is fixed at any value z. 

At each longitudinal coordinate z„ the radiation field can be expanded in a 
series of modes of the passive (empty) waveguide. Although we have not 
calculated explicitly the coefficients in this expansion, these are independent 
of z, as z = Zi has been fixed. This leads to a paradox. In fact, one can 
propagate each waveguide mode till another coordinate Zf. However, each 
mode of the passive waveguide has its own phase velocity. Therefore, one 
will not recover Eq. (87) at position z/, because relative phases of modes 
have changed, at that position. The paradox is solved noting that the field 
in Eq. (87) represents a result for an active (loaded) waveguide. In other 
words, to obtain back Eq. (87), one should add to the propagated field, the 
additional field radiated by the electron beam between the points z, and 
Zf, which is nothing but the edge radiation contribution from the setup in 
Fig. (19) with A = Zj and B = Z/. In other words, in order to recover the 
self-reproducing field structure, we should add to the right hand side of the 
propagation equation 
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where G is given in Eq. (81). Note that Eq. (88) generalizes Eq. (14). It follows 
that 
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/ dz'[d'fl{?fj\zf-z')\,__^exp 

Note that the second term on the right hand side of Eq. (89) is nothing but 

the edge radiation contribution from z, to Zf (for z, = -L/2 and Zf = L/2 one 
may compare with Eq. (78)). Eq. (89) can be rewritten as 



2/ 

£ER(L/2,r) = ^ j dz' [d'^G"^{rf,f,Zf-z')l,^^exp 

Zi 

= £y L/2) - ^ J dfE^,{-L/2, r^)G;{K r\ L) , (90) 

where we defined Esi(-L/2) = -Es_i.ep(-L/2)andEs2(L/2) = Es_i.ep(L/2). Hav- 
ing introduced this notation we are in the position to see that the treatment 
of edge radiation in terms of virtual sources has a straightforward gener- 
alization for the case of a waveguide. In fact, Eq. (90) means that the edge 
radiation field can be interpreted as the sum of a given field distribution 
E"^^L/2) at z = L/2 plus the field obtained propagating another field dis- 
tribution, E"^{-L/2,f), from z = -L/2 up to z = L/2 according to Eq. (88). 

We thus straightforwardly interpret £si,s2 as the generalized virtual sources 
at z = -L/2 and z = L/2 respectively. Note that in the limit for p — > oo, 
the self-reproducing field solution in Eq. (87) is related to the virtual source 
fields in Eq. (39) consistently with what has just been stated. 

Also note that Eq. (90) and Eq. (84) present the same field, and must coincide. 
This hints to the fact that the sum in Eq. (84) can be calculated analytically. 
Since the use of virtual sources provides conceptual insight and should 
facilitate the design of ER setups, it is instructive to explicitly demonstrate 
the equivalence of Eq. (90) and Eq. (84). This is demonstrated in Appendix 
A. 



7.3 Analysis of results (exemplifications) 

We now analyze our main results, Eq. (84) and Eq. (85). To this purpose, it 
is convenient to re-write these equations in normalized units, so that: 
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(92) 



where we defined = C^L = Vq^/(2Q), we reminded Q = p^/(AL), and 
normalized quantities E, f and z are defined as in unbounded space. 

Note that, with reference to Fig. 19, we are assuming z < 1/2, because for 
z > 1 /2 we are already in unbounded space. Using the Green's function 
we can formally propagate the field at z = 1/2, i.e. Eq. (84), up to positions 
z < 1 /2. Note that, since Eq. (84) is presented as a sum of empty-waveguide 
modes, this propagation just modifies the exponential (phase) factor, so that 
instead of exp[-zQ72] we find exp[-zC^z] in Eq. (91). Then, Eq. (91) assumes 
the meaning of a virtual field distribution in analogy with the unbounded- 
space case. 

Also note that the sinc( ) function in the expression for is a direct conse- 
quence of our model of the setup of sharp-edges at positions z = -1/2 and 
z - 1/2 respectively. In fact, the sinc(-) function is the Fourier transform of a 
rectangular function with respect to C^ + (p/ 2, modelling the sharp switch on 
and switch off of the bunch harmonics. The presence of high frequency com- 
ponents in the rectangular function implies the presence of contributions 
with high values of k. In its turn, each contribution with large value of k in 
Eq. (91) can be interpreted as a superposition of plane waves propagating 
at angles Vqa:/ VQ in units of the diffraction angle yfXfh. Thus, higher values 
of k correspond to the introduction of high spatial frequency components in 
the field. 

However, we should account for the fact that our theory applies with a finite 
accuracy related to the value of the parameter 5. When using this approxi- 
mation we neglect contributions to the field with an accuracy related with 
b. This means that, in the analysis of our results, it does not make sense to 
consider high spatial frequency contributions due to abrupt switching of 
the bending magnet fields on a scale shorter than ^/fFa, because these are 
outside of the accuracy of the sharp-edge approximation. We may then in- 
troduce a spatial frequency filter in our expression for the field by replacing 



Eq. (78) with 
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where the function S(z') introduces some smoothing of the rectangular pro- 
file on a scale "V^R^A. In normalized units Eq. (93) reads 
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£«=47i r dz'S{z')[d'i^G''Jlf',z-z')] exp 



(94) 



where derivatives are now calculated with respect to normalized quantities. 

We model S(z') as a constant function along the straight section length with 
exponentially decaying edge^ on a typical (normalized) distance A: 



exp [-(z' + 1/2)2/(2A2)] for z' < -1/2 
S(z') = < 1 for - 1/2 < z' < 1/2 

exp [-(z' - 1/2)V(2A2)] for z' > 1/2 

It follows that Eq. (92) should be replaced by 
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where !F{S(z'), (C[' + (p/2)} is the Fourier transform of the function S with 
respect to (Cj[''^ + (p/2). The introduction of the function S(z') introduces a 
suppression of higher frequency components of !F{S}, and corresponds to 
a suppression of higher spatial frequencies in the field distribution, i.e. to 
a smoother field distribution. For a fixed value of Q and z we qualitatively 
expect finite and smooth field distribution for finite values of A. On the 
contrary, in the sharp-edge limit for A — > 0, the field must diverge for some 
value of f, as a result of the already discussed divergence of the integrated 
flux and of the fact that the radius of the vacuum chamber is finite. Note 
that the level of high spatial frequencies depends in a complicated way on 
Q, because the perfectly metallic waveguide effectively acts as a mirror. The 



^ Formally, such model is in contrast with the previous statement that for z > 1/2 
we are in unbounded space. However, we introduce A to account for a finite 
accuracy of the sharp-edge approximation. Thus, it does not make sense to specify 
where the unbounded space begins within a distance A. 
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Fig. 21 . Comparison between hard-edge case and high-frequency filtering. Function 
S(z') in the hard-edge case with A = (a) and in the case when filtering is applied 
with A = 0.1 (b). Their Fourier transforms TiS} with respect to ^ = C^^ + C>Q are 
plotted in (c) for the hard edge case and in (d) when filtering is present. 




Fig. 22. Intensity profiles of the virtual source (z = 0) at large values of Q = R^/(ALa,) 
for different values of A in the limit for ^ <sc 1. These plots are obtained from Eq. 
(91) and Eq. (96), while the asymptotic limit for A «: 1 is found with the help of Eq. 
(31). 

value of A should be actually chosen to cut off high spatial frequencies that 
are outside the region of applicability of the sharp-edge approximation, i.e. 
A ~ 6. Thus, the correct value of A depends, case by case, on the type of 
setup considered. The Fourier transform of the function S(z') in Eq. (95) to 
be inserted into Eq. (96) reads 
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where the imaginary error function erfi(-) is defined as 
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Fig. 21 presents a comparison between functions S and !F{S} for A = 0.1 and 
A = 0. A filtering effect can be clearly seen, suppressing higher frequency 
components of 'F{S}. 

Eq. (91) and Eq. (92) completely describe ER in the presence of a circular 
waveguide. Position 2 = corresponds to the single virtual source descrip- 
tion discussed before. In fact, as we have seen, the concept of laser-like ra- 
diation beam can be naturally extended to the case a waveguide is present. 
As said before one must carefully select the number of modes used for 
computation, according to A: » VQ. 

We further consider, as an example, the limit for (/) «; 1. In Fig. 22 we show a 
comparison between the analytic expression for the intensity distribution I 
at the virtual source position, and numerical expressions obtained through 
Eq. (91) for different values of A at Q = 1000 and Q = 100. Differences 
between results for A = 0.1 and Eq. (31) should be taken as exemplification 
of the filtering process. A three-dimensional view of the intensity profile 
of the virtual source for Q = 5 and A = 0.1 is shown in Fig. 23, to be 
compared with Fig. 9 for the unbounded-space limit. The number of modes 
used in all plots are in all cases much larger than VQ. Furthermore, we 
verified that results do not change by changing the number of modes used 
in the computation, provided that condition k^ax » is fulfilled. We also 
checked that the same results obtained with the help of Eq. (91) can also 
be recovered by introducing apodization directly in Eq. (19). This confirms 
that numerical evaluation of Eq. (91) is sound. 

The behavior of the intensity profile at the virtual source when waveguide 
effects are important is shown in Fig. 24 for different values of Q at A = 0.01 
A = 0.1. Further on, in Fig. 25 we plot the intensity distribution I = |Ep as a 
function of the reduced radius f for the same values of A and Q at z = 1/2. 
All plots are evaluated with the help of Eq. (91). It should be remarked that 
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Fig. 23. Three-dimensional plot of the intensity profile of the virtual source (z = 0) 
at Q = K^I{XLy,) = 5 for A = 0.1 in the limit for «: 1. 

the wavefront at z = 1 /2 is not plane (as at z = 0) and the phase distribution 
carries important information when one needs to propagate the field from 
z = 1/2 to any other point. The derived analytical algorithm in Eq. (91) 
enables one to calculate all these informations. 

Comparison of the cases for A = 0.01 and A = 0.1 clearly shows the ef- 
fect of different edge dimensions, corresponding to filtering of high spatial 
frequencies. It should be reminded that when A = the angle-integrated 
spectral flux of ER is divergent in the unbounded space, as well as in a 
waveguide. Such divergency was motivated as an effect of the sharp-edge 
approximation, and depends logarithmically on A. 

To sum up, knowledge of the field at z = 1/2 fulfill our goal for this Sec- 
tion. Expression for the field at z = 1/2 can be interpreted as a virtual-field 
distribution. From this viewpoint, the presence of a finite parameter A cor- 
responds to the presence of an angular filter, related to the violation of ER 
expression for angles comparable with the bending magnet radiation angle. 
The virtual-field distribution at z = 1/2 should be further propagated either 
in the unbounded space or through any other kind of vacuum chamber up 
to the experimental station. In other words, the knowledge of the field at 
the straight section exit constitutes the starting point for further propaga- 
tion. Typically, a user-beamline vacuum chamber will be connected with the 
electron vacuum chamber (that follows the bending magnet after position 
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A=0.1 



A=0.01 




Fig. 24. Intensity profiles of the virtual source (z = 0) at different values of 
Q = R2/(ALjy) for ^ 1 at two values A = 0.1 and A = 0.01. 



z = 1/2). The user-beamline vacuum chamber has a typical characteristic 
size which is larger than the electron vacuum chamber. It follows that, after 
position z = 1/2, the field can usually be propagated as in unbounded space. 
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Fig. 25. Intensity profiles at z = 1/2 at different values of Q = R^/{XL„,) for ^ «c 1 at 
two values A = 0.1 and A = 0.01. 

8 Coherent edge radiation 

The physics of radiation processes from electron beams critically depends 
on the ratio of the bunch length to the wavelength of the emitted radiation. 
We can thus define two opposite asymptotes for the electron bunch length, 
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corresponding to very different characteristics of radiation. 

In the first limit, when the electron bunch is long compared to the radiation 
wavelength, one deals with the conventional case of spontaneous radiation. 
In this regime the electron phases are not correlated, i.e. distances between 
electrons are randomly distributed on a characteristic scale much longer 
than the radiation wavelength. Incoherent superposition of each electron 
contribution to the radiation field yields low-power emission, proportional 
to the charge in the bunch, with the usual incoherent phase-noise statistics. 

In the second limit, when the electron bunch is short compared to the ra- 
diation wavelength, one deals with the case of coherent radiation. Electron 
phases are correlated, in the sense that random variations of distances be- 
tween electrons are distributed on a characteristic scale shorter than the ra- 
diation wavelength. The bunch essentially behaves as a single point-charge, 
which coherently emits high-power radiation. In this case, the radiation 
power scales quadratically with the total electric charge in the bunch, rather 
than linearly as in the other limit. It is because of this quadratic effect that 
the radiation power can be greatly enhanced, because the bunch population 
is usually of order 10^'^ particles. 

Coherent radiation is usually not emitted in electron storage rings, because 
the bunch length is of order of a centimeter. However, in linear accelerators, 
much shorter electron bunches of the order of lOOfim can be produced in 
magnetic bunch compressors. In this Section we will discuss the physics of 
coherent ER distinguishing two cases, when coherent emission is due to an 
overall temporal bunch profile shorter than the wavelength or to short struc- 
tures (microbunches of order of the radiation wavelengtih.) superimposed to 
the temporal profile of the bunch. 

When it comes to the use of ER as a diagnostic tool, the availability of a much 
larger number of photons constitutes an important advantage of coherent 
ER compared to the incoherent case. By detecting coherent ER, character- 
istics of the electron bunch such as their three-dimensional distribution 
and divergence can be measured. Moreover, coherent radiation presents the 
unique feature that electron-beam microbunching can be investigated too. 

8.1 Review of known methods 

There are two equivalent approaches to describe coherent SR in general: a 
Lagrangian approach and a Eulerian approach. 

In the Lagrangian approach particles are labelled with a given index, and 
the motion of individual charges is followed through space. One follows the 
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evolution of each particle as a function of energy deviation, angular direc- 
tion, position and arrival time at a given reference (longitudinal) position. 
Knowing the evolution of each particle, individual contributions to the field 
are separately calculated and summed up. 

In the Eulerian approach instead, one begins defining charge density p{r, z, t) 
and current density /(r, z, t), i.e. particular field quantities, as a function of 
coordinates and time. One is no more interested in knowing the evolution 
of each single particle but only to observe charges passing through a given 
longitudinal position, in order to obtain information about p and /. Once 
these two fields are determined, one can solve Maxwell's equations for the 
electric field, with p and / as macroscopic electromagnetic sources. 

8.1.1 Lagrangian approach 

We consider an electron bunch moving along the setup specified in Fig. 1(a). 

Given a reference point, e.g. z - 0, along the trajectory one may specify 
the phase-space coordinates of a given electron with the help of six data: 

position of the electron, /, angular deflection rj, energy 6y, and arrival tim.e, 
T, all relative to the position of a reference electron with position ro(0), 
transverse velocity Uxo(O), nominal energy Sq and arrival time f = 0, all 
given at position z = 0. Using the index k to identify a certain particle 
one finds the temporal Fourier transform of the radiation pulse at a given 
frequency cu as a sum of many independent contributions: 



where Ne is the total number of particles in the bunch, and E(i) is the 
single particle field contribution from a given particle identified by vari- 



ables /fc, 6y]^, Tfc. Note that we separated these variables from the space- 
frequency coordinates where the field is calculated with a semicolon. 

Eq. (99) forms the basis for calculating the angular spectral flux of radiation. 
This is essentially done by taking the square modulus of the field in Eq. (99) 
and averaging over an ensemble of bunches: 




(99) 



k=l 




(100) 
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Let us show that in the particular case when no focusing elements are 

present, E(i) in the far zone depends on through a phase only. To this 

purpose, it is sufficient to perform the substitution ro(z') — > J'o(z') + I in 

-> 

Eq. (13), which leads straightforwardly to the phase -cod ■ l/c. Thus, in the 
particular case when there are no focusing elements and electrons have 
random arrival times (phases) and random offsets, but no energy spread, 
nor deflections, Eq. (100) yields 
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having recognized that if ft is the longitudinal bunch profile and // the offset 
distribution at the reference position z = 0, one has: 



00 

(expiicoTk)) = J dzk ft{Tk) exp{icoTk) = ft{co) 
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and 
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In the rest of this Section we will mainly concentrate on the second term in 
Eq. (101), scaling as ~ NJNe - 1), i.e. the coherent term. Although Eq. (101) is 



known (see e.g. i36|.l37|. l38(]). we choose to shortly report it, because it is still 
the object of revisions in recently published articles |39l, l40l, 14111 , where the 
angular spectral flux of spontaneous Transition Radiation is found to obey 
a modified angular distribution depending on the transverse distribution of 
electrons (see e.g. Eq. (25) of jilll , dealing with the case of an infinitely large 
mirror). In the next Section we will demonstrate that setups in Fig. 1(a) and 
Fig. 1 (c) can be treated with the help of the same mathematical machinery, i.e. 
a Transition Radiation setup can be formally dealt with as if it were an edge 
radiation setup. While an exponential cutoff is to be expected in the coherent 
angular spectral flux, results in | j39|. 40|. 41] indicating a dependence of the 
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incoherent angular spectral flux on the transverse distribution of electrons 
is counterintuitive. In fact, as we have seen, the spectrum is calculated by 
first summing single electron contributions to the field and, second, by 
averaging the square modulus of the field over an ensemble of bunches. As 
we have discussed above, a given offset modifies the far-field contribution 
from a smgle electron by a phase factor only. Our conclusion, in contrast 



with |39l,l4d,l41|1 , is that the final incoherent spectrum cannot depend on the 
transverse distribution of electrons. The incorrectness of [39,4fl/i41J is due to 
a calculus mistake, which can be found e.g. in Eq. (19) of [4l']: the average of 
the square modulus of a quantity, in fact, is not equal to the square modulus 
of the average. 

Let us go back to the most general situation described in Eq. (99). One can 
represent the sum in Eq. (99) in terms of an integral over a continuous phase- 
space distribution of electrons specified at a given position z. Let us assume 
for simplicity that this reference position is z = 0, i.e. in the middle of the 
straight section. We will indicate the phase-space distribution at z = with 

f{l, if, 5y, t; 0). Eq. (99) can then be substituted with an average of the single 
particle field over /, yielding 



E (r, z, co) - j dldfjd{6y)dT f (l, if, by, E(i) {l, ff, by, t; r, z, a;) . 

(104) 

Note that substitution of Eq. (99) with Eq. (104) rules out the possibility of re- 
covering the incoherent term in Eq. (101). In fact, the continuous phase-space 
distribution / already includes an average over an ensemble of bunches. In 
other words, Eq. (104) has the meaning of a mean field, whose square mod- 
ulus is essentially the second term in Eq. (101), i.e. the coherent term. 



8.1.2 Eulerian approach 

Up to now we have been using a microscopic approach, where the single- 
electron field is averaged over the six-dimensional phase-space distribution 
of electrons. It is possible to obtain the same results using a macroscopic 
viewpoint, where solution of Maxwell's equations is found with given 
macroscopic charge density p(r, z, t) and current density /(r, z, t). In this case, 
one can integrate Maxwell's equations with the help of the usual Green's 
function. 

We begin by writing a general expression for the temporal Fourier trans- 
form of charge and current density in terms of the phase-space distribution 
at a given position down the beamline z. Such phase space-distribution. 
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f{K^±/c,dy,t,z) can always be written as /(/ + fo{z),ff + v_L_o{z)/c,6y,T + 
s{z)/v,z) where, as we have seen before, the subscript "0" referring to the 
constrained motion of an electron with r = 0, Vj^ = 0, 6y = and t = at 
z - 0. The meaning of this writing is that here we consider the evolution of 
/ along z as the result of two combined effects. On the one hand we have a 
z-dependent shift of the /-distribution as a whole due to the reference mo- 
tion ro{z), v±o{z) and s(z), while on the other hand we have a more generic 
evolution in z due to self-interactions and external fields. The distribution / 

is obviously of the form f{l, t], by, t, z). It follows that the Fourier transform 
of /, /, can be written as 



f{l,ri,6y,co,z) - I dtexp[ia)t]f{l,ri,6y,T,z) 

oo 




f{l,r],6y,co,z) , 



(105) 



where the usefulness in the definition of / is that / is a slowly function of z 
on the scale of A. One has therefore 



p(r, z, co) = ^ ^^^^ dldi]d{by)f{l, ff, by, co, z)b{l + ro(z) - r ) 



i-e)N, 



J difdiby) f{r- ro(z), rf, by, co, z) exp 
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and 



j{r, z, co) = ^ J" dldifdiby) (crj + Vj_{z)) f{l, ff, by, co, z)b{l + ro(z) - r) 



c 

(-g)N, 
c 



J drfd{by) (cff + v±{z)) f{r- ro(z), ff, by, co, z) exp 

icos{z) 



= j{r-ro{z),z,co) exp 



ia)s{z) 

V 

(107) 



where we introduced, together with / and with the usual field envelope £ = 

£exp[-za;z/c], also quantities p = pexp[-za;s(z)/u] and; = ;exp[-za;s(z)/y]. 
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Note that p = p{f- ro(2)/ z, co) (the same applies to ;). One writes the paraxial 
equation for the field as 



— -r-iE-4n exp < ico 

c dz) \ \ V ci 



X 



ICO- 



-Yi(r- ro{z), z, co) - V^p (r - ro(z), z, co) 



L c- 



(108) 



in analogy with Eq. (9) for a single particle. We solve Eq. (108) with the 
help of the Green's function in Eq. (10), similarly as in Eq. (11), but with 
the important difference that, now, we cannot carry out the integration over 
transverse coordinates. We therefore obtain 



z 

'=471 J dz' J dr'^^%r -rQ{z'),z',co)G[r-r',z-z') 

—oo 

+p{f - ro(z'), z', cv)V\G (f- f, z - z')) exp jzo; |^ - ^ 
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Explicit substitution of Eq. (10) in Eq. (109) finally yields 



E(z,r)- I dz I dr' exp < ico — -+ | 

^ ^ c J j z - z' ^ I 2c(z - z') j 



dz 



lyl{z)c 
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c z-z' 
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Lagrangian and Eulerian approaches are equivalent. In simple cases, e.g. 
when trajectories of electrons are a simple function like a straight motion and 
Eq. (13) can be calculated analytically, the Lagrangian method looks more 
transparent. However, in the most general case, one needs to account for 
the presence of focusing elements and self -interactions. The evolution of the 
six-dimensional phase space becomes complicated and numerical methods 
should be used. Naturally, the Lagrangian approach can be applied when 
the evolution of the six-dimensional phase-space is known via simulations. 
Then, Eq. (13) has to be calculated numerically for each (macro)particle. The 
Eulerian approach can be more advantageous. Computer codes can calculate 
charge and current density distribution as a function of z accounting for 
focusing elements and self -interactions, and the macroscopic method may 
be more practical, because one integrates Eq. (110) once and for all. One does 
not need to have any explicit knowledge of the angle distribution, or of the 



73 



energy spread distribution to calculate the field. Once charge and current 
density are known, one only needs to (numerically) integrate Maxwell's 
equations. This is because we are dealing with a macroscopic system, where 
only macroscopic current and charge density are important. In other words, 
the behavior of Maxwell's equations is influenced by projections of the 6D 
phase-space on the real space (p and /), but not by the 6D phase-space itself. 



8.2 Near Field Coherent Edge Radiation 



Let us consider the usual setup in Fig. 1(a) and assume, till the end of the 
Section, that there are no focusing elements nor self -interactions along the 
straight section. By this, we assume that electrons move along the straight 
line between two bending magnets. We will adopt the Lagrangian approach 
to characterize coherent ER emission due to overall temporal bunch profile, 
without considering, for now, microbunching. 

An algorithm for describing radiation emission is straightforward. One 
begins with the expression for the field generated in the far zone by a single 
electron occupying a certain position in phase-space at a fixed reference 
longitudinal position z. Then, one integrates the single-particle expression 
over the phase-space distribution, thus getting the equivalent of Eq. (104) in 
the far zone. Finally, when the far-zone field (averaged over a given electron 
distribution) is known, the field in the near-zone can be calculated as has 
been done before, by back-propagating the far-zone expression. 

In the case of edge radiation, the far-zone expression for the field of a single- 
particle having phase-space coordinates (/, rj, 6y, t) at z = is given by: 



£(1) (l, f], by, t; 6, z, co^ - exp \i(x)%\ exp 

x(0-rf)sinc|^ 



icod^z 
2c 

d-i] 



exp 



icoO • I 
c 

26y 



y2 y3 



. (Ill) 



Note that Eq. (Ill) is a straightforward generalization of Eq. (20), obtained 
by performing a rigid rotation of the reference system and a translation of 
the observation point without changing the observation plane of interest, 
which remains perpendicular to the z axis of the initial reference system. 

This explains why the combination d — rj does not enter the phase in Eq. 
(111). This solves the problem for the field in the case when trajectories 
are straight lines between two bending magnets without focusing elements, 
and self-interactions are negligible. 
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The field in the far zone should be subsequently averaged over /. This 
is accomplished by substituting Eq. (Ill) into Eq. (104). Simplifications 
arise when the phase-space distribution can be separated in product of 
different factors depending only on offset, deflection, energy spread and 
arrival time for some choice of the reference position z. Such separability is 
not always granted. However, the choice of some particular reference point 
may help simplifying calculations. For example, the distribution of offsets 
and deflections factorizes when the reference position in z is fixed at the 
minimial value of the betatron function. In the following we will consider 
two different study cases of interest. 



I. Case / = fi(l)6{rf)6{6y /'y)fj{T). This is one of the simplest situations, 
when both angular distribution and energy distribution can be neglected, 
typical of well-collimated, monochromatic high-quality beams produced in 

conventional accelerators. Simpler subcases for /; = 6{l) or fj = 6(t) (or 
both) can be easily inferred from this more general model. 

As we discussed before, the presence of an offset modifies the single-particle 
field by a phase factor only. In case of ER, Eq. (104) simplifies to 



£(z,0,a;) =N,f (a;,0)E(i) 
icoeNeFio), 6)1 



exp 



icod^z 
2c 




y2 



(112) 



where 



F(co,d)^ J dtf,{z)exp[icoT] J dlfi{l)exp -y^-?] • (113) 

Thus, in this case, we obtain a simple result amounting to a multiplication 
of the structure factor F defined in Eq. (113) by the single-particle field. The 
structure factor F is the product of the Fourier transform of fj (the bunch 
longitudinal profile), fj, and the Fourier transform of the bunch transverse 
profile. 

To sum up, if one knows results for a single electron, like those calculated 
above for ER, TUR, or any other kind of radiation, and one needs to calculate 
radiation from an electron bunch, all is needed to do is to multiply the result 
for a single electron by f ; results in the near zone are obtained remembering 
that 6 = r/z. Note that in the case of a single electron with zero offset and 
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_ —T—T 

arrival time = 0, = 1 and //(/) = so that F = 1 and one recovers Eq. 
(13). 

In the present case I, the form factor F does not change during the motion of 
electrons along the z axis. This is in agreement with our previous assumption 
about energy spread and angular distribution of the bunch. In fact, as said 
before, we assume that /t(t) remains unchanged during the motion along 
the beamline. Concerning the transverse structure, we are assuming that 
the transverse size and shape of the bunch do not vary either. 

With the help of Eq. (24) and Eq. (112) we obtain the angular spectral flux 



= |F(a;,0)| 



2 cNlz^ 

47l2 



(114) 



Note that Eq. (114) is just the second (coherent) term in Eq. (101). We should 
stress that the near-zone problem for coherent ER can be solved, similarly 
as in the single particle case, starting from the far-zone data. Therefore, it 
can be solved using Eq. (112). 

Consider first the case ^ 1. The problem of calculation of the virtual 
source distribution in the middle of the straight section can be basically 
reduced to a Fourier transform of the far-field distribution. Therefore, the 
specification of virtual sources also simplifies to a convolution between the 
single particle virtual source and a spatial Fourier Transform of the form 
factor F. Remembering Eq. (16), and with the help of Eq. (20), Eq. (30) and 
Eq. (112) one explicitly obtains 



£ (0, r) = N,/;(a;)£(i)(0, f) * Mr) 

Aicoe 



= NeMco) j 



dl 



smc 



CO 



2m 



cL 
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(115) 



where £(i)(0, indicates the single-particle field distribution at the virtual 
source in Eq. (30), and the symbol " * " denotes the convolution product as 
defined in the second line of Eq. (115). Note that the virtual source in Eq. 
(115) can be imaged by a lens when the object plane corresponds to position 
z = 0. 



Similar remarks hold for the case of a generic value of (/). In fact, as is 
well-known from Eq. (37), the expression for the single-particle field in 
the far-zone can be split in the sum of two spherical-wave contributions. 
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and both terms are subsequeritly multiplied by the same form factor. Thus, 
analogously to Eq. (115), one obtains 



£ (+L/2, r) = Ne/;(a;)£(i)(±L/2, r) * //(r) , (116) 
-» 

where now E(i)(±L/2, r) indicate the single-particle field distribution at the 
virtual source positions ±L/2 in Eq. (39). 



11. Case / = fi,ri{LTf)H^y/y)f'z{'^)- In this situation we still neglect the 
energy-spread distribution, but we account for a given distribution in an- 
gles in addition to the case treated in I. In other words, the bunch has a finite 
emittance. Factorization of /; and fjj at the reference point (z = 0) indicates 
that the betatron functions are minimal at that position. In practice, angular 
distribution is important for the case depicted in Fig. 1(d), i.e. for a laser- 
plasma accelerator device. In this case, it is more natural to assume that the 
minimal betatron function is located at z = -L/2, rather than at z = 0. Then, 
the product fi{l)frj{f]) must be consistently substituted by a non-factorized 

distribution //,,,(/, rf). We will consider a Gaussian electron bunch with rms 
size CT;. Additionally, we introduce a Gaussian angular distribution with 
rms and assume, for simplicity, that these rms quantities apply both in 
the horizontal and in the vertical plane at z = -L/2. As a result, one has 
the following distribution at z = (which may be alternatively rewritten in 
terms of Twiss parameters): 



{I - Lrj/2f t]' 
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(117) 



With this in mind, Eq. (104) simplifies to 



E{z,d,(jo^ = 



icoeNeLf^ico) 
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Generalization of this study case to the situation with a finite energy spread, 
still relevant e.g. for laser-plasma accelerator devices can be performed 
straightforwardly. 
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8.3 Microbunching 



We will now study the particular case of an electron bunch modulated in 
density at a given wavelength. In this case, coherent emission is not due to 
the overall bunch profile, but to the presence of a modulation superimposed 
to such profile. 

In some cases modulation of the bunch can be set up on purpose, while in 
others it can be a detrimental effect. 



Typical examples of parasitic microbunching mechanisms at PEL facili- 
ties are driven by Coherent Synchrotron Radiation (CSR) in the bending 



magnet of a bunch compressor chicane |42l, |43|, |44|1 or Longitudinal Space 
Charge force (LSC) [45]. Initial shot-noise density modulation at a given 
frequency induces energy modulation at the same frequency due to CSR 
or LSC impedances. Such energy modulation is subsequently transformed 
back into density modulation through a dispersive section, eventually lead- 
ing to a beam instability. The dispersive section is typically the magnetic 
compression chicane itself. In this case, microbunching is usually an un- 
wanted phenomenon, which tends to spoil the high-quality electron beam 
needed for the PEL process. 

However, in other cases, a modulation at optical wavelengths can be printed 
on purpose onto the electron bunch with the help of a (quantum) laser and a 
modulator undulator followed by a dispersive section. Optically modulated 
electron bunches can be used for diagnostic purposes and in pump-probe 



techniques as well ||46l,l24|1 



If the bunch is modulated at a given wavelength A,„, one has /t(t) = 
/to(t)[1 -I- flo cos(cT/Am)], where /^o is a given longitudinal profile and Uq 
is the modulation level. 



Let us consider a typical Gaussian model for the beam, specified by 



(119) 



Ox being the rms beam transverse dimension and by 



/.o(t) 



exp 



24. 



(120) 



Here we are assuming that the wavefront of the density distribution of the 
modulation is perpendicular to the velocity of electrons, and is uniform in 
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the transverse direction as Uq is not a function of r. The latter assumption is 
close to reality, because the waist of the seed laser is much larger than the 
size of the electron bunch in the modulator undulator, due to the extreme 
high-quality of the electron bunch in XFELs. 

Near the modulation frequency cOm = c/Am, a bunch with modulated Gaus- 
sian temporal profile and rms duration Oj yields 



M<^) = y I exp 



\2 

a>m) 



+ exp 



(121) 



where /^(w) is, as before, the Fourier transform of the temporal profile of the 
bunch. In deriving Eq. (121) we used an adiabatic approximation that can 
be taken advantage of in practical situations involving XFELs, where the 
optical modulation wavelength is much shorter than the bunch length, i.e. 
cojIAm » 1. In this case we have another large parameter in the problem, 
which considerably simplifies the treatment. In fact, since we are interested 
in coherent emission around optical wavelengths, in calculating Eq. (121) we 
neglected the Fourier transform of the first term of /, (t), i.e. /to(t^)- Obviously, 
this would not be possible for frequencies in the range II Ot- From Eq. (121) 
follows that radiation is exponentially suppressed for frequencies outside 
the bandwidth Ijoj centered at the modulation frequency co^- 

The exponential suppression in Eq. (121) has its analogous in the transverse 
direction too. With the help of Eq. (119) one finds 



F [ci), 0) = /t(£l') exp 



(122) 



Note that Q is usually normalized to a characteristic angle that depends 
on the system considered. In our ER case, such angle is ^KjLj, with — 
min [L, y^A]. With this in mind, the exponential suppression in Eq. (122) can 
be written in normalized units as 



Fa; 



exp 



(123) 



where we defined 



(124) 
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The parameter N is analogous to a Fresnel number in diffraction theory, and 
is the only (dimensionless) transverse parameter related to the radiation 
emission in the model defined by Eq. (119) and Eq. (120). When N ^ 1 
the exponential factor in Eq. (123) returns unity, and one gets back the 
single-electron radiation profile, whereas the flux is possibly modified by 
the presence of When N > 1 the intensity profile as a function of angles 
is exponentially suppressed. 

Integration of Eq. (114) in dco and in dO, gives the total energy radiated. Since 
we are dealing with an electron bunch modulated at optical wavelengths 
we are interested in narrow-bandwidth emission. 

Within the above-mentioned adiabatic approximation we can consider the 
wavelength in E(i) appearing in Eq. (114) fixed, i.e. A = Am. Then, inte- 
gration in dco can be easily performed using the fact that dco |/T(ct')|^ = 
yjnay {^Oj), leading to 



do. 



exp 



(125) 



Integration over angles can be performed analytically in particular limiting 
cases. 



In the following we will discuss the setup in Fig. 1(a). However, our result 
can be easily extended to the case of Fig. 1(b) and Fig. 1(c). Let us first 
consider the case cp <^1. One has 
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oo 
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Introducing dimensionless units allows one to write the total number of 
photons emitted at cOm as 



oo 

^vh = — -T-sin^ ^ exp = ' ° In 1 + — ,127) 



where a = e^/{hc) = 1/137 is the fine structure constant and Nm = OjCOm is 
the number of modulation wavelengths included in the bunch. It should be 
noted that Eq. (127) is logarithmically divergent for N — > 0. It is therefore 
interesting to discuss the region of applicability of Eq. (127) in the limit 
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for small values of N. Such asymptote can usually be exploited when deal- 
ing with XFEL setups, because of the extremely high-quality of the electron 

bunch (small emittance). The expression used for E(i) assumes that the sharp- 
edge approximation holds, i.e. 6 ^ 1. In this case, the formation length of 
the bending magnet at the end of the straight section Lfb is much smaller 
than the straight section length L. The transverse spot-size of the radiation 
associated with bending magnet emission, and is of order ALfb is then the 
smallest characteristic transverse size of the problem. If the electron bunch 
size becomes of order AL^, our electrodynamic description cannot distin- 
guish anymore between a bunch with finite transverse size and a point. In 
this case one should substitute oi with ALp. This amounts to the substitution 
N — > 6. Thus in the limit N «; 6 we prescribe the substitution of N with 6 
for estimations. 

Finally, since cotI Km » 1, it makes sense to introduce an expression for the 
instantaneous power as a function of the instantaneous (peak) current /(t) 
and modulation aoiz), if the seed laser duration is not too long compared to 
the bunch duration. To this purpose, we consider a stepped-profile model 
for the bunch /to(t) - 1/T for -T/2 < t < T/2 and zero elsewhere. It follows 
that 
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Then, in the limit for T » A^/c one has dco |/t(£l')|^ - alnT, and 
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where dP/ dQ, is the angular distribution of power. If now the current I = I(t) 
and the modulation level a = a{T) are slowly varying functions of time on 
the scale A^/c, one obtains for cf) 
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(130) 



Eq. (130) can be used in order to study the general case of an electron bunch 
with arbitrary gradient profile and amplitude of modulation. 

Let us now turn to the case when » 1 for the single-edge radiation 
{y^A <^ dz ^ L and r <^L/y). Introducing apodization of the detector on a 
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characteristic angular scale Oa and using Eq. (42) we obtair0 



oo 
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r ,,2 c2^2 







(131) 



Note that apodization (through Oa) and transverse bunch length (through 
^m/oi) play the same role from a mathematical viewpoint. Also note that, 
in order to use Eq. (42), we must require that we are within its limits of 
applicability, i.e. rfi • iTiin(crfl, Am/oi) L/y (always with y'^Am « ^2 ^ L). 

Taking advantage of dimensionless units one can write the total number of 
photons emitted at oj,„ as 



where A = A/{2Lol) is analogous to the N parameter in Eq. (124) and ac- 
counts for apodization effects, and T{a,z) = t"~ exp{—t)dt is the incom- 
plete gamma function. Discussion of asymptotic cases for different values 
of N and A is possible on the basis of Eq. (132). 



9 Extraction of edge radiation by a mirror 

9.1 Backward Transition and Diffraction radiation 

A setup where edge radiation from an upstream bending magnet is ex- 
tracted by a mirror is shown in Fig. 1(c). In the particular case depicted in 
that figure, radiation is reflected by a hole mirror and sent to a far-infrared 
spectrometer. This example refers to a method | j47l ] where coherent radiation 
is used to monitor the bunch length at XFEL setup, but similar setups can 
also be used to collect incoherent radiation in the optical range. The size and 
shape of the mirror may vary. A hole may be present in the mirror or not, 

^ In particular, Eq. (131) applies to the case when coherent TR due to microbunching 
is outcoupled by a mirror and postmirror optics can be modelled with a Gaussian 
pupil profile. In this case, z - L/2 indicates the distance from the mirror. 



oo 




[-1 + [l + (A + N)(/)] r[0, (A + N)(/)] exp[(A + N)(p\] ,(132) 
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Fig. 26. Geometry of backward transition radiation (BTR) (a), backward diffraction 
radiation (BDR) (b). 

and the geometry of the experimental arrangement depends on the partic- 
ular setup considered. Anyway, in its basic lines, a standard edge-radiation 
diagnostic setup consists, similarly as in Fig. 1(c), of a mirror positioned at 
some distance after a bend and rotated by an angle n /4 in the direction of the 
electron bunch to allow extraction of radiation through a vacuum window. 
Such setup can be used for longitudinal bunch-length measurements in the 
mm and sub-mm range as well as for transverse electron bunch diagnostics 



in the optical wavelength range |47l,l48l,l49l] . 



Consider the straight section after the bend. If this is long enough so that 
(/) » 1 (and 6 ■ (p « 1, i.e. A » Ac), radiation collected by the mirror is 



the usual Transition Radiation (TR) first noted by Ginzburg and Frank [5C1( ] 
and considered in a number of publications written during the last half 
century. In recent years there has been a good deal of interest in TR, because 
this form of radiation has found useful applications, e.g. in diagnostic for 
ultrarelativistic beams [51.,.5Z,.53.,,54,,55,] . These papers describe an electron 
crossing an interface between two media with different dielectric constants 
and they refer, in particular, to the case of the boundary between vacuum 
and an ideal conductor. As a consequence of the crossing, time-varying 
currents are induced at the boundary. These currents are responsible for 
TR (see Fig. 26 left). The metallic mirror, that is treated as the source of 
TR, is usually modelled with the help of a Physical Optics approach. This 
is a well-known high-frequency approximation technique, often used in 
the analysis of electromagnetic waves scattered from large metallic objects. 
Surface current entering as the source term in the propagation equations 
of the scattered field are calculated by assuming that the magnetic field 
induced on the surface of the object can be characterized using Geometrical 
Optics, i.e. assuming that the surface is locally replaced, at each point, by its 
tangent plane. 

As mentioned before, one typically assumes that the Ginzburg-Frank for- 
mula can be used, i.e. the field distribution is proportional to Ki[a)r/ (yc)]. In 
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terms of parameters introduced in this paper, this approximation has sense 
only when » 1. However, even for this case, one should always specify 
the transverse region of applicability of the usual TR formulas, which is 
typically neglected (as we will see, this region is for r <^ (pAy = L/y, with 
(/) » 1). Moreover, mainly due to an increased electron energy, today setups 
often meet, in practice, condition (/) < 1 even in the optical wavelength 
range. In this case one should account for the whole setup when discussing 
the reflected radiation, and not only for the metallic mirror. In other words, 
as we will see, the problem of extraction of edge radiation reduces to the 
well-studied TR problem only asymptotically, i.e. for » 1. In this paper 
we will still keep on talking about TR or, with some abuse of language, even 
when usual TR formulas do not apply anymore ^ 1/ or r > L/y). 

When dealing with the above-discussed setup, one usually talks about Back- 
ward Transition Radiation (BTR). BTR is widely used for different purposes, 
because it allows for low background radiation levels. We will consider a 
whole range of setups for BTR, without treating their geometry in detail, 
and we will give an algorithm to deal with these systems in all generality. In 
fact, in all cases, the main problem is in the specification of the electric field 
distribution at some position where a mirror is present (the hole mirror in 
the example of Fig. 1(c)). The only element to be accounted for, aside for the 
mirror itself, is an upstream dipole magnet. Note that the mirror is tilted 
as in Fig. 1(c), i.e. it is not perpendicular with respect to the optical axis. In 
principle, one may account for this tilting, but in practice the projection of 
the mirror on the optical axis is negligible compared with the distance to 
the dipole, and one may consider the screen as not tilted. 

Note that Diffraction Radiation (DR) appears when charged particles move 
in the vicinity of a medium (e.g. a conductive mirror) at some given impact 
parameter ^, and has been recently suggested as a possible tool for non- 
invasive bunch diagnostics. DR is obviously related to TR and, in particular, 
it can be treated starting from the knowledge of the electric field distribution 
at the mirror, and subsequently applying physical optics techniques. BDR 
in the optical wavelength range has been measured and applied for trans- 
verse electron beam diagnostics (see e.g. [48, 49]). A schematic comparison 
between a BDR and a BTR emission is shown in Fig. 26. 

It should be stressed that specification of the field at the target position 
must be considered as the first step to the solution of a more complicated 
problem, i.e. the characterization of the field at a detector position. Such 
first step is considered separately, because the field at the target plane is 
independent of the type of target and detector. Once the field at the target 
position is known, the full problem can be solved with the help of Physical 
Optics techniques, where one should account for diffraction effects due to 
the particular shape of the mirror. 
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In the following we will focus on the problem of field characterization at 
the target position. A virtual source method can be applied. Our algorithm 
consists of the following steps: (i) propagate the field from the upstream 
virtual source to the mirror, (ii) add the field from the downstream virtual 
source. 

Note that, in contrast to the study case of a straight section between two 
bends considered before, numerical methods involving direct integration of 
Maxwell's equations (e.g. with the help of SRW) down the optical axis up to 
the mirror will fail, because electromagnetic sources must be propagated up 
to the observation plane, and integration of the Green's function up to that 
position yields a singularity (independently of the value of 6). In our virtual 
source approach we also deal with a singularity, but this is isolated in a 
single term, and expressed by means of an analytical function (the Bessel Ki 
function in Eq. (39)), allowing straightforward presentation of results. The 



same approach has been proposed in ||56[1 to discuss the setup in Fig. 1(c). 



Let us focus in detail on our generic BTR setup. Following the algorithm 
described above, we model the system with the help of two virtual sources at 
z = +L/2, i.e. at the end of the bend and at the mirror. As said before, Eq. (39) 
can be used to describe the virtual sources. Application of the propagation 
formula Eq. (14) allows to calculate the field at any distance z in free-space 
for the source located at z = -L/2. In order to simplify the presentation of 
the electric field we take advantage of polar coordinates. We propagate the 
upstream virtual source (located at z = -L/2) up to the mirror position at 
z = L/2, and we add the field distribution for the downstream source at 
z = L/2. 

At z = L/2 (i.e. at the mirror position) we obtain 
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which is valid for any value of (/). Exactly as in the case treated in Section 5, 
one may deal with two asymptotes of the theory for » 1 and <^1. 

In the case (/) » 1 the contribution from the upstream source in Eq. (133) 
can be calculated analytically. In fact, the Ki function under the integral sign 
is responsible for an exponential suppression of the integrand in the limit 
for r' » yA. As a result, we can integrate up to values r' < yA without 
changing the integration result. This means that the phase factor in the 
integrand will be at most of order Ijcp « 1, and can be neglected. The result 
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of the integration then corresponds to the expression for a single edge in 
the far-zone, Eq. (41), the only difference being that the single-edge source 
is located at z = -L/2. Altogether, for (p » 1, Eq. (133) simplifies to: 
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Eq. (134) is valid for r <^L i/A/R, i.e. for angles (measured from the upstream 
edge at z = -L/2) smaller than the angle of SR from the bend. When r « yA, 
the Ki term is obviously dominant with respect to the second (polynomial) 
term in r. When r > yA, the Ki term drops exponentially, while the second 
term is zero at r = 0, grows polynomially and reaches its maximum at r ~ 
L/y. Since the polynomial growth is slower than the exponential decrease of 
the term, one will have a region yA < r ^ L/y where the field is strongly 
suppressed with respect to the maximum at r ~ L/y. In this region, the Ki 
term is still a good approximation to the total field. As a result, the first 
term in Eq. (134) can be taken as a good approximation for the total field 
when r ^ L/y (always assuming cp » 1). In this asymptotic case, results 
from our analysis coincide with usual results from TR theory. Note that the 
question about how many photons will be collected strongly depends on 
the acceptance of the optical system. The first term in Eq. (134) is singular 
at r = 0, so that the larger the acceptance, the nearest to the singularity we 
collect photons, and the largest the number of photons is. 

The opposite asymptote for (p -^1 can be described by applying the propa- 
gation formula, Eq. (14), to Eq. (30), like it was done for Eq. (32). When z is 
fixed at the mirror position, i.e. at z = L/2 one obtains |57l] 
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The quadratic phase factor in Eq. (135) describes a spherical wavefront 
centered on the optical axis at the upstream edge (the initial bend). Note that 
this contribution scales as 1 /r and is singular on-axis. The screen positioned 
at the downstream end of the setup will detect electric field given by the 
spherical wavefront in Eq. (135) at any value y/AR^ <c L « y^A from the 
upstream magnet and r <^L -s/A/R. 



Following the specification of the field (in amplitude and phase) at the tar- 
get position, one faces the problem of propagation of BTR in free-space 
and through optical elements, which has been extensively discussed in 
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|5^|55(]. As mentioned above, this reduces to a standard Physical 



Optics problem, which is interesting within the field of Optical Engineering, 
i.e. from the viewpoint of the design of post-mirror collector optics. In these 
references, discussions are based on the assumption that a field distribution 
Ki[cor/{yc)], i.e. the Ginzburg-Frank formula, can be used. In this case, we 
deal with a laser-like beam, whose waist is located on the target (i.e. the 
target exhibits a plane wavefront). The transverse size of the waist is of or- 
der yA » A, so that Fourier Optics can be applied to solve the propagation 
problem. The Rayleigh length of this laser-like beam is about Ri - y^K, 
and the Fresnel diffraction-zone is located at distances from the mirror of 
order Rl- When the distance from the mirror becomes much larger than 
one has the Fraunhoffer diffraction-zone. If some linear optical element is 
located at a given position after the mirror, one can use available wavefront 
propagation/physical optics codes (e.g. SRW 125], ZEMAX [58J, GLAD [59], 
PHASE [60]) for calculations. In other words, the propagation problem is 
well-defined, and does not include any novel aspect from a scientific view- 
point. Only, one needs to put attention to special terminology and notations 
used in literature, which may lead to some misunderstanding. For example, 
the Rayleigh length R^ ~ y'^K is named "the formation length". Since the 
waist is always given by Wq ^ -y/L^, it follows that Rl ~ vo^^jK = Lf. How- 
ever, it should be noted that, conceptually, the formation length can only be 
discussed in connection with the particle trajectories, and not in free-space. 
Similarly, the Fresnel diffraction-zone is usually named as the "pre-wave 
zone", while the Fraunhoffer diffraction-zone is called "wave-zone".From 
this viewpoint, articles dealing with pre-wave-zone effects actually treat a 
routine Optical Engineering problem. 

It should be remarked that the problem of specifying the field at a mirror 
position is fundamentally different from what we have done in the previous 
Sections. Before we considered a system composed by a straight section 
between two magnets, and radiation was detected downstream of the last 
magnet. That setup could be treated entirely within a SR-theory approach, 
and computer codes like SRW could be used to cross-check our results. Here 
the situation is different. In fact, Eq. (11) must now include integration up 
to a position z where sources are present 
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-L/2 



with z = L/2, and therefore includes a singularity. In analogy with Eq. 
(90) one then obtains Eq. (133), where the singularity is isolated in the first 
field-term, which is entangled with sources, whereas the second term is 
disentangled. 
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Note that in the mathematical limit for p — > oo, Eq. (90) reduces to Eq. 
(133). Similarly, Eq. (40) reduces to Eq. (133) at z = L/2. However, although 
formally similar, Eq. (90) and Eq. (40) have different physical meaning com- 
pared to Eq. (133). In this Section we are characterizing the field along the 
straight section. This gives a different status e.g. to the singularity in Eq. (133) 
for r = and z = L/2 compared to the singularity in Eq. (40), also located 
at r = and z = L/2. In fact, the singularity in Eq. (40) (which reduces to 
Eq. (32) in the case (p ^ 1) is related with the sharp-edge approximation: 
6 has a finite value, and going beyond the accuracy of the sharp-edge ap- 
proximation (i.e. modelling the edges and accounting for their contribution 
to the radiation) the integrated flux will not diverge, nor the reconstructed 
field from the far zone data will. In contrast to this, the singularity in the 
Ki term in Eq. (133) at r = and z = L/2 is due to the fact that the plane at 
z = L/2 intercepts the electron trajectory at that point, and is fundamentally 
related with our model of the electron beam as a filament without transverse 
dimensions. 



9.2 Impact on setups at XFEL facilities 

During the last decade, the electron beam energy in linac-based facilities 
increased up to around 10 GeV and novel effects, which were considered 
negligibly small at lower energies, became important. In fact, at high ener- 
gies, condition L » y^A becomes unpractical even for optical wavelengths, 
and the effect of the presence of an upstream bending magnet on TR emis- 
sions must be considered. Moreover, when measuring the long-wavelength 
coherent BTR, i.e. when the radiation wavelength is comparable to, or longer 
than the electron bunch length, the influence of the upstream bending mag- 
net could be significant even at lower energies. These facts led to several 
misconceptions that can be found till recently. 

An example of application for diagnostic purposes in XFEL setups is given 
in (^^. In general, monitoring methods based on coherent radiation take 
advantage of the sensitivity of the coherent spectrum on the bunch form 
factor, which is related to the bunch length. Coherent radiation from a given 
setup is thus collected with the help of mirrors and its energy measured, for 
example with the help of pyroelectric detectors. In order to analyze results, 
one has to characterize the field distribution at the collecting mirror position. 

The setup proposed in ||47tl uses the last bend in each magnetic compression 
chicane of LCLS as depicted in Fig. 27. Radiation is collected by a mirror 
with a hole (that has to be large enough to let the electron beam through) and 
imaged to the diagnostic station. The hole mirror is positioned after a short 
straight section, a few tens of centimeters long. Detailed parameters can be 
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y = 500 ' imaging 

R = 2.2 m 

lull ~ 50 cm-' 
mirror is located 16 cm 

after the last dipole 



diaenostic 



Fig. 27. Layout of the LCLS bunch length monitor, located immediately after the 
last dipole of the first magnetic bunch compressor. 

found in ||47|]. In the case of the BCl chicane, for example, the length of the 
straight section between magnet and mirror is L = 0.16 m, the wavelength 
of interest is of ordeJZJA ~ 1 mm, the electron energy is about 250 MeV, and 
the radius of the bend is about R - 2.2 m. 



In |l47t] one can read: "the critical wavelength of the synchrotron radiation is 
orders of magnitude smaller than the wavelength range considered here for 
the bunch length monitor and the synchrotron radiation can be neglected 
compared to edge radiation". Authors of |0] stress the fact that A » Ac. Such 
condition is absolutely necessary to satisfy the applicability of ER theory. 
However, as we have seen, it is not sufficient, as the extra-requirement 6 <^c 1 

must be satisfied too. If we calculate the parameter 6 = y/R^A/L for numbers 
given above we obtain 6 - 0.6. We conclude that synchrotron radiation from 
the bend cannot be neglected: on the contrary, its effect will be as important 
as those from the straight section. The sharp-edge approximation fails here. 

Note that there are two contributions to the field at the mirror position. 
The first is the contribution due to bending magnet radiation, while the 
second is the contribution due to the straight section. The fact that 6 ^ 
0.6 indicates that both contributions are of the same order of magnitude. 
It follows that it is incorrect to neglect the bending magnet contribution. 
Assuming that calculation of the straight section contribution is performed 
correctly, this would lead to a quantitatively incorrect numerical result from 
the viewpoint of practical applications. However, authors of ||47|l calculate 
the straight section contribution incorrectly as well. Let us consider the (p 
parameter. The reader may check that in this case the formation length 
y^A - 50 m is much longer than the length of the straight section L = 0.16 m, 
so that (p <^ 1. Therefore, Eq. (135) should be used to calculate the straight- 
section contribution at the mirror position. Authors of ||47l] proceed instead 



Since we are talking about coherent radiation, the reduced characteristic wave- 
length A is related here with the longitudinal size of the bunch Orms/ i-e. A ^ 
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Fig. 28. ODR setup at the APS facility. In this experiment, ODR at visible wavelength 
is emitted for impact parameter ^ between 1 mm and 2 mm, values close to the 
Lorentz factor y times the reduced observation wavelength. 

in the following way. First, they recognize that for ER "the Coulomb field 
is the radiation source due to the sudden change in acceleration". Second, 
they propagate a source from the edge of the magnet up to the mirror 
"by expressing the Coulomb field as a sum of Gauss-Laguerre modes". 
This corresponds to the second term of Eq. (133). In other words, the field 
distribution at the mirror position is obtained by propagating a K\ source 
located at the upstream bending magnet. It follows that the result in ||47l] is 
qualitatively incorrect. 



To sum up, in 1471], SR from the bending magnet is incorrectly neglected, and 
in the calculation of ER the downstream source, the first term of Eq. (133), 
is missed, qualitatively changing the result. The correct procedure to find 
the field distribution at the mirror for this parameter values (5 ~ 1, (/) « 1) 
is to use Eq. (135) to calculate the straight-section contribution at the mirror 
position, and to sum up the SR contribution from the bending magnet. 



We now turn to analyze |48l, |49(], where a setup similar to that in Fig. 28 
is studied from an experimental viewpoint. An electron bunch travels a 
straight section of length L - 5 m downstream of a bending magnet before 
intercepting an offset mirror at the ODR station. The wavelength considered 
here is A - 800 nm, and the electron energy is 7 GeV. In this case y'^K - 30 
m, i.e. (/) - 0.17. Authors of |48, use a downstream virtual source ~ 
]<Ci[^/(yA)] to characterize the field at the BDR screen, effectively describing 
an ER setup by means of usual Gtnzburg-Frank expression [50] for TR. 
However, here 6^1 and «: 1. As we have seen before, in this parameter 
range the correct procedure to evaluate the field distribution at the mirror is 
to use Eq. (135), which is different from the Ginzburg-Frank expression. In 
other words, in this parameter-range, the Gtnzburg-Frank formula cannot 
be applied anymore. 



Note that in |48|1 is reported that "the ODR peak signal intensity depen- 
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dence is consistent with exp[-2rf/(yA)], where d is the distance from beam 
center to the screen edge, y is the Lorentz factor". However, in this region 
of parameters the asymptotic behavior in Eq. (135), does not exhibit any 
exponential cutoff, nor dependence on yA. Therefore, the interpretation of 
experimental result in [48] cannot be accepted. 

This example and that in [47] help to understand the importance of edge- 
radiation theory in relation with future XFEL facilities. They also demon- 
strate the usefulness of similarity techniques in the efficient planning of 
experimens. 



10 Conclusions 

In this article we showed how the theory of laser beams can be used to char- 
acterize radiation field associated with any Edge Radiation (ER) setup. In 
fact, in the space-frequency domain, ER beams can be described in terms of 
laser-like beams, with large transverse dimensions compared to the wave- 
length. Similarly to usual laser beams, ER beams were shown to exhibit a 
virtual "waist" with a plane wavefront. The field distribution of ER across 
the waist turned out to be strictly related to the inverse Fourier transform of 
the angular field-distribution in the far-zone. As a result, standard Fourier 
Optics techniques could be taken advantage of, and the field could be prop- 
agated to characterize ER beams at any position down the beamline. In 
particular, we reconstructed the near-field distribution from the knowledge 
of the far-field ER pattern. This could be accomplished by (i) describing 
the far-field pattern with known analytical formulas, (ii) finding the virtual 
source(s) and (iii) propagating the virtual source distribution in the near 
(and far) zone. 

After a qualitative discussion in Section 2, we applied our techniques to a 
typical setup constituted by a straight section between bends in Sections 3, 4 
and 5. These Sections constitute the first comprehensive treatment of ER, in 
the sense that we consistently used similarity techniques for the first time, 
allowing discussion and physical understanding of many asymptotes of the 
parameter space together with their region of applicability. In particular, the 
main parameters of the theory are found to be 6, the ratio between the bends 
formation length and the straight section length, and (p, the ratio between 
the length of the straight section and the maximal formation length of ER. 
Note that introduction of the parameter 6 allowed us to define for the first 
time "how sharp" the edges are, and to specify the region of applicability 
of ER theory. A classification of regions of observation of interest, which is 
regarded by us as a novel result, is presented in Section 5.3 with the help of 
dimensionless parameters. 



91 



In Section 6 we applied our treatment to deal with a Transition Undulator 
Radiation (TUR) setup. As before, we relied on virtual source expressions 
derived from the far-field pattern. These virtual sources were propagated 
in free-space in the near zone, thus providing for the first time an exact 
analytical characterization of TUR in the near zone. 

In the following Section 7 we presented the first exhaustive theory of ER 
within a waveguide. The analysis of the problem was performed by intro- 
ducing a tensor Green's function technique. This complicates the mathemat- 
ical structure of equations that depends, contrarily to the unbounded-space 
case, on the waveguide geometry. We outlined a solution for a homoge- 
neous waveguide with arbitrary cross-section, further specializing it to the 
case of a circular waveguide. The electric field was found as a superposition 
of waveguide modes, and was studied for different values of parameters. 
The main parameter involved in the problem (other than (p, 6 and z) was 
found to be a waveguide parameter Q, related with the strength of guiding 
effects. This can be interpreted (at (/) «: 1) as the ratio between the waveg- 
uide radius and the radiation diffraction size, and is a purely geometrical 
parameter. 

In Section 8 we presented considerations on the coherent emission of ER. The 
availability of a large number of photons constitutes an important advantage 
of coherent ER compared to the incoherent case when this type of radiation 
is used as a diagnostic tool, where characteristics of the electron bunch can 
be measured, and electron-beam microbunching can be investigated too. 

Finally, in Section 9 we discussed the problem of extraction of ER by a 
mirror, stressing differences and analogies that this kind of setup presents, 
compared to those discussed in the previous Sections. Particular attention 
was given to diagnostics setups for XFEL facilities. 

As a final remark, it should be noted that in our work we consistently 
exploited both theoretical and numerical results from simulation. These 
approaches are complementary, and we first took advantage of such com- 
plementarity. On the one hand, there are situations when existing codes 
cannot be applied, e.g. the case of a waveguide (Section 7) and the char- 
acterization of the field distribution at a mirror position (Section 9), where 
part of the field is entangled with the sources. In this case, analytical results 
can be used for practical calculations. On the other hand, computer codes 
can easily account for finite bending magnet edge length (finite value of 6) 
in a particular set of problem parameters. From this viewpoint, our theory 
can be used to prepare, based on similarity techniques, particular sets of 
problem parameters to be used as input for computer codes, which subse- 
quently return universal plots presenting the accuracy of ER theory in terms 
of dimensionless parameters. 
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Appendix: Equivalence of Eq. (90) and Eq. (84) 



We want to show that Eq. (90), together with Eq. (87), is equivalent to Eq. 
(84). We can demonstrate this fact by calculating the sum in Eq. (84), and 
presenting it in terms of an integration. 

We begin expanding the sinc(X) = sin(X) / X function in terms of exponential 
functions. Moreover, we subtract from X in the denominator of sin(X)/X 
an infinitesimally small quantity ie, where e > 0, to be dropped at the 
end of calculations. In other words, we substitute 1/X with 1/(X - ie) = 

Jr»00 
P dE, exp[-zX<^ - e<^]. We therefore obtain 
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Eq. (137) consists of two field terms, corresponding to the two terms inside 
the {...} parenthesis. Let us consider the first term and demonstrate that 
such term is in fact a generalization of the virtual source at z = L/2 (i.e. 
Eq. (39) for Zg = L/2) accounting for the presence of the waveguide. Using 
the fact that /i(X) = l'Q{iX)li, the ""' symbol indicating derivative with 
respect to the argument, and remembering the Wronskian relation 1/X = 
/o(X)Ki(X) + /i(X)Ko(X), we have 
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Using Eq. (138) and remembering /i(X) = Ii{iX)/i we write the first term in 
Eq. (137) as 
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We now consider the integral 
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where Oi is an integration path on the complex r/ plane going from y' - zoo 
to y' + ioo, y' > Q being a positive number larger than the real part of all 
singularities of the integrand, and closed by a semicircle of infinite radius 
in the left half plane. It can be shown that the integrand is a single valued 
function of rj. It follows that u can be calculated as a sum of residues of first 
order poles located at r] = r]jc = -iALv'^J{4p'^). We obtain 
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. r;=-!cLv2 /(4a;p2) 



that can be immediately substituted in Eq. (139). The integral along the 
semicircle in Eq. (140) can be dropped using Jordan's lemma. Therefore 
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Eq. (142) is easily seen to be the Laplace transform of an inverse Laplace 
transform. Letting e — > 0, it follows straightforwardly that 
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Eq. (143) is a generalization for the virtual source located at z = L/2 that 
accounts for the presence of a waveguide (compare with Eq. (87)). 

Let us now consider the second term in Eq. (137), which can also be written 
as 



exp 



(144) 



Since 



I 



+Voj2(Vofc)Xi — 

_ cop \cy) 



-—Ji{vok)K2\ — 

cy \cyj 



(145) 



and 



/ 



[0)0 IcopX {cop\ 
^^^j/iW+Vo./.(^jHvo.) 



(146) 



one can write Eq. (144) as 
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(147) 

because the second and the third terms on the right hand side of Eq. (145) 
cancel with Eq. (146) multiplied by Kq (a)p/{cy))/Io (a)p/{cy)). Remembering 
that an expression for the field propagation from some initial point z, to 
some final point Z/ is given by Eq. (88), one finds that Eq. (147) is obtained 
by propagating the field distribution 
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up to Zf - L/2. This can be seen with the help of Eq. (81), selecting m - 0. 

Now, Eq. (148) is just a generalized expression, valid in the presence of a 
waveguide, for a virtual source at z = -L/2, which is described similarly to 
Eq. (143), the difference being only its location position. 

It follows that the second term in Eq. (137) is obtained propagating the 
virtual source at z = -L/2 up to position z = L/2, that is the second term of 
Eq. (90). 
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